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Abstract This article provides an overview of various 
notions of shape spaces, including the space of parame- 
trized and unparametrized curves, the space of immer- 
sions, the diffeomorphism group and the space of Ric- 
mannian metrics. We discuss the Riemannian metrics 
that can be defined thereon, and what is known about 
the properties of these metrics. We put particular em- 
phasis on the induced geodesic distance, the geodesic 
equation and its well-posedness, geodesic and metric 
completeness and properties of the curvature. 
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1 Introduction 

The variability of a certain class of shapes is of in- 
terest in various fields of applied mathematics and it 
is of particular importance in the field of computa- 
tional anatomy. In mathematics and computer vision, 
shapes have been represented in many different ways: 
point clouds, surfaces or images are only some exam- 
ples. These shape spaces arc inherently non-linear. As 
an example, consider the shape space of all surfaces of 
a certain dimension and genus. Therefore it is difficult 
to do statistics. One way to overcome this difficulty 
is to introduce a Riemannian structure on the space of 
shapes. This may allow to locally linearize the space and 
develop statistics based on geodesic methods. Another 
advantage of the Riemannian setting for shape analysis 
is its intuitive notion of similarity. Namely, two shapes 
that differ only by a small deformation are regarded as 
similar to each other. 

In this article we will concentrate on shape spaces 
of surfaces and we will give an overview of the different 
Riemannian structures, that have been considered on 
this space. 

1.1 Spaces of interest 

We fix a compact manifold M without boundary of di- 
mension d — 1. In this paper a shape is a submanifold 
of M d that is diffeomorphic to M and we denote by 
B z (M,R d ) and B e (M,R d ) the spaces of all immersed 
and embedded submanifolds. 
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Diff(M) s- lmm(M,R d ) ^ Diff c (R d ) 




Vol+(M)/Diff(M) Met(M)/Diff(A/) Met(R d ) 

Fig. 1 An overview of the relations between the spaces discussed in this article. Where Vol_|_(M) is the space of all positive 
definite volume densities, Diff(M) the diffeomorphism group, Met(M) the space of all Riemannian metrics, Imm(M, l* 1 ) the 
space of immersed surfaces and B;(M, R d ) the shape space of unparametrized surfaces; g denotes the Euclidean metric on R d . 



One way to represent Bi(M,Mr) is to start with 
the space Imm(Af, R ) of immersions upon which the 
diffeomorphism group Diff(M) acts from the right via 



A way to define Riemannian metrics on the space 
of parametrized submanifolds is via the left action of 
Diff c (M d ) on Emb(Af,R d ), 



Imm(M,R d )xDiff(Af) 3 (q,<p) ^qotpe lmm(M,R d ). Diff c (R d ) x Emb(M, R d ) 3 (tp, q) H> <p o q <= Emb(M,] 



The space Imm(Af, Mr) is the space of parametrized 
type M submanifolds of R d and the action of Diff(AI) 
represents reparametrizations. The immersions q and 
q o ip have the same image in R d and thus one can iden- 
tify Bi(M,R d ) with the quotient 

B,{M,R d ) = Imm(M,R d )/Diff(M) . 

The space Bi(M, R d ) is not a manifold, but an orb- 



ifold with isolated singular points; see Sect. |3.3| To re- 
move these we will work with the slightly smaller space 
Imm/(Af, R d ) of free immersions and denote the quo- 
tient by 

Bij(M,R d ) S Imm/(M,R d )/Diff(Af) . 

Similarly one obtains the manifold B e (M, R d ) as the 
quotient of the space Emb(Af, R d ) of embedded sub- 
manifolds with the group Diff(Af). 

The spaces Imm(Af, R d ) and Emb(M, R d ) are some- 
times called pre-shape spaces, and Diff(Af) is the re- 
parametrization group. Their quotients Bi(M,R d ) and 
B e (M, R d ) are called shape spaces. 

A Riemannian metric on Imm(Af, Mr) that is in- 
variant under the action of Diff(Af) induces a metric 
on B i j(M,R d ), such that the projection 

7T : Imm/(Af,R d ) -> B itf (M,R d ), q i-> q(M) , 

is a Riemannian submersion. This will be our method 
of choice to study almost local and Sobolev-type met- 
rics on lmm(M,R d ) and B i j(M ) M d ) in Sect. [5] and |6| 
These are classes of metrics, that are defined via quan- 
tities measured directly on the submanifold. We might 
call them inner metrics to distinguish them from outer 
metrics, which we will describe next. This is however 
more a conceptual distinction rather than a rigorously 
mathematical one. 



Given a right-invariant metric on Diff c (W. d ), the left ac- 
tion induces a metric on Emb(M, R d ), such that for 
each embedding qo £ Emb(M , R d ) the map 

TT qo : Diff c (M d ) -> Emb(Af,R d ), ip.^ipoq , 

is a Riemannian submersion onto the image. This con- 
struction formalizes the idea of measuring the cost of 
deforming a shape as the minimal cost of deforming the 
ambient space, i.e., 



GT b {h,h) 



inf 

X o q=h 



^.Diff 



(X,X). 



Here h g T q Emb(M, R d ) is an infinitesimal deforma- 
tion of q and the length squared mb (h,h), which mea- 
sures its cost, is given as the infimum of Gf d (X,X), 
that is the cost of deforming the ambient space. The 
infimum is taken over all X € X c (M d ) infinitesimal de- 
formations of R d , that equal h when restricted to q. 
This motivates the name outer metrics, since they are 
defined in terms of deformations of the ambient space. 

The natural space to define these metrics is the 
space of embeddings instead of immersions, because not 
all orbits of the Diff c (IR' i ) action on Imm(Af , R rf ) are 
open. Defining a Riemannian metric on B e {M, R d ) is 
now a two step process 



-Emb(M, R d )—^-B e (M, R d ) 



Diff, 



One starts with a right-invariant Riemannian metric on 
Diff c (R d ), which descends to a metric on Emb(Af,R d ). 
This metric is invariant under the reparametrization 
group and thus descends to a metric on B e (M, R d ). 
These metrics are studied in Sect. [8] 

Riemannian metrics on the diffeomorphism groups 
DiS c (M. d ) and Diff(Af) are of interest not only, because 
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these groups act as the deformation group of the ambi- 
ent space and the reparametrization group respectively. 
They are related to the configuration spaces for hydro- 
dynamics and various PDEs arising in physics can be in- 
terpreted as geodesic equations on the diffeomorphism 
group. While a geodesic on Diff c (IR d ) is a curve <p(i) of 
diffcomorphisms, its right-logarithmic derivative u{t) = 
d t (fi(t) o </?(i) _1 is a curve of vector fields. If the metric 
on Diff c (K d ) is given as Gid(u, v) — J Rd (Lu, v) dx with 
a differential operator L, then the geodesic equation can 
be written in terms of u(t) as 

<9fTO + (u ■ V)m + to div u + Du T .m — 0, m — Lu . 

PDEs that are special cases of this equation include the 
Camassa-Holm equation, the Hunter-Saxton equation 
and others. See Sect. for details. 

So far we encoded shape through the way it lies in 
the ambient space; i.e., either as a map q : M — > M. d or 
as its image q(M). One can also look at how the map q 
deforms the model space M. Denote by g the Euclidean 
metric on K d and consider the pull-back map 

Imm(M,R d ) -> Met(M), q h-> q*g , (1) 

where Met(M) is the space of all Riemannian metrics 
on M and q*g denotes the pull-back of g to a metric on 
M. Depending on the dimension of M one can expect 
to capture more or less information about shape with 
this map. Elements of Met(M) with dim(M) = d - 1 
are symmetric, positive definite tensor fields of type ( 2 ) 
and thus have d ^ d ~ 1 ^ 1 components. Immersions on the 
other hand are maps from M into M. d and thus have 
d components. For d — 3, the case of surfaces in IR 3 , 
the number of components coincide, while for d > 3 we 
have > d. Thus we would expect the pull-back 

map to capture most aspects of shape. The pull-back is 
equivariant with respect to Diff(M) and thus we have 
the commutative diagram 



Imm/(M, K d ) 



Bij{M,R d ) 



q^q g 



Met(M) 



Met(M)/Diff(M) 



The space in the lower right corner is not far away from 
Met(M)/Diff (M), where Diff (M) denotes the con- 
nected component of the identity. This space, known as 
super space, is used in general relativity. Little is known 
about the properties of the pull-back map ([!]) , but as a 
first step it is of interest to consider Riemannian met- 
rics on the space Met(M). This is being done in Sect. 
11 with special emphasis on the L 2 - or Ebin-metric. 



1.2 Questions discussed 

After having explained the spaces, that will play the 
main roles in the paper and the relationships between 
them, what are the questions that we will be asking? 
The questions are motivated by applications to com- 
paring shapes. 

After equipping the space with a Riemannian met- 
ric, the simplest way to compare shapes is by looking 
at the matrix of pairwise distances, measured with the 
induced geodesic distance function. Thus an important 
question is, whether the geodesic distance function is 
point-separating, that is whether for two distinct shapes 
Co 7^ C\ we have d(C^,C\) > 0. In finite dimensions 
the answer to this question is always "yes" . Even more, 
a standard result of Riemannian geometry states that 
the topology induced by the geodesic distance coin- 
cides with the manifold topology. In infinite dimensions, 
when the manifold is equipped with a weak Rieman- 
nian metric, this is in general not true any more. The 
topology induced by the geodesic distance will also be 
weaker than the manifold topology. We will therefore 
survey what is known about the geodesic distance and 
the topology it induces. 

The path realizing the distance between two shapes 
is, if it exists, a geodesic. So it is natural to look at 
the geodesic equation on the manifold. In finite dimen- 
sions the geodesic equation is an ODE, the initial value 
problem for geodesies can be solved, at least for short 
times, and the solution depends smoothly on the initial 
data. The manifolds of interest in this paper are natu- 
rally modelled as Frechet manifolds and in coordinates 
the geodesic equation is usually a PDE or tJTJE. Thus 
even the short-time solvability of the initial-value prob- 
lem is a non-trivial question. For some of the metrics, 
in particular for the class of almost local metrics, it is 
still open. For the diffeomorphism group the geodesic 
equations for various metrics are of interest in their own 
right. To reflect this we will discuss in Sect. [7] first the 
geodesic equations before proceeding with the proper- 
ties of the geodesic distance. 

It is desirable for applications that the Riemannian 
metric possesses some completeness properties. It can 
be either in form of geodesic completeness, i.e., that 
geodesies are extendable for all time, or metric com- 
pleteness with respect to the induced geodesic distance. 
Since we are considering only weak Riemannian metrics 
on spaces of smooth shapes, we cannot expect the space 
to be metrically complete, but in some examples it is 
possible to at least describe the metric completion of 
shape space. 

In order to perform statistics on shape space one can 
choose or compute a mean shape and linearize the space 
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around this shape via the Riemannian exponential map 
and normal coordinates. The curvature tensor contains 
information about the accuracy of this approximation. 
In general computing the curvature leads to long for- 
mulas that are hard to interpret, but in some cases 
curvature admits a simple expression. We collect the 
examples, where more is known about the curvature, 
usually the sectional curvature, than just the formula 
for it. 

To summarize, we will deal with the following four 
properties of Riemannian metrics on shape spaces: 

— Geodesic distance 

— Geodesic equation and existence of geodesies 

— Geodesic and metric completeness 

— Properties of the curvature 

Riemannian geometry on shape spaces is currently 
an active area of research. Therefore this paper is less 
an encyclopedic treatment of the subject but rather re- 
sembles an interim report highlighting what is known 
and more importantly, what is not. 



1.3 Topics not discussed 

There are many topics that lie outside the scope of this 
paper, among which are the following. 

Changes in topology. In certain applications it may 
be of interest to consider deformations of a shape that 
allow for the development of holes or allow the shape 
to split into several components. In this paper we fix 
the model manifold M and only consider submanifolds 
of R d diffeomorphic to M. Thus by definition all de- 
formations are topology-preserving. See 20J40 122 for 
topologically robust approaches to shape matching. 

Non-geodesic distances. Many interesting distances 
can be defined on shape spaces, that are not induced 
by an underlying Riemannian metric; see for example 
As we are looking at shape spaces through 
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the lens of Riemannian geometry, these metrics will nec- 
essarily be left out of focus. 

Subgroups of the diffeomorphism groups. The Rie- 
mannian geometry of the diffeomorphism group and its 
subgroups, especially the group of volume-preserving 
diffeomorphisms, has been studied extensively; see for 
example 113 . It plays an important role in hydrody- 



namics, being the configuration space for incompress- 



ible fluid flow 42 . While the full diffeomorphism group 



itself is indispensable for shape analysis, its subgroups 
have not been used much in this context. 

Utmost generality. We did not strive to state the 
results in the most general setting. It is possible to con- 
sider shapes of higher codimension inside Mr or curved 



like space curves or curves that lie on a sphere. It would 
also make the presentation more difficult to read. 

Numerical methods. Since shape space is infinite- 
dimensional, computing the exponential map, the geo- 
desic between two shapes or the geodesic distance are 
numerically non-trivial tasks. While we present some 
examples, we do not attempt to provide a comprehen- 
sive survey of the numerical methods that have been 
employed in the context of shape spaces. Finding sta- 
ble, robust and fast numerical methods and proving 
their convergence is an area of active research for most 
of the metrics and spaces discussed in this paper. See 
[4p7|38 ,55,105:107 for various approaches to discretiz- 
ing shape space. 

2 Preliminaries 

2.1 Notation 

In this section we will introduce the basic notation that 
we will use throughout this article. On M. d we consider 
the Euclidean metric, which we will denote by g or •). 
Unless stated otherwise we will assume that the param- 
eter space M is a compact, orientable manifold without 
boundary of dimension dim(Af) = d — 1. Riemannian 
metrics on M are usually denoted by g. Tensor fields 
like g and its variations h are identified with their as- 
sociated mappings TM —> T*M. For a metric g this 
yields the musical isomorphisms 



b : TM -> T*M and 



T*M -> TM. 



Immersions from M to R d - i.e., smooth mappings with 
everywhere injective derivatives - are denoted by q and 
the corresponding unit normal field of an (orientable) 
immersion q is denoted by n q . For every immersion q : 
M ~R d we consider the induced pull-back metric g — 
q*g on M given by 

q*g(X,Y)=g(T q X,T q Y), 

for vector fields X, Y £ X(M). We will denote the in- 
duced volume form of the metric g — q*g as vol(g). In 
coordinates (it, U) it is given by 



vol(g) = ^Jdet(g (d^fyqj) Idu 1 A • • • A du d-1 | . 

Using the volume form we can calculate the total vol- 
ume Vol g = J M vo\(q*g) of the immersion q. 

The covariant derivative of a metric g will be de- 
noted by V 9 and we will consider the induced Bochner- 
Laplacian A 9 , which is defined for all vector fields X 6 
X(M) via 



ambient spaces; see [ll]. This would include examples A 9 X = Tr(g- 1 V 2 X) . 
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Note that in R d the usual Laplacian A is the negative 
of the Bochner-Laplacian of the Euclidean metric, i.e., 
A3 = -A. 

Furthermore, we will need the second fundamental 
form s q (X, Y) — g {v q x ^Tq.Y,n^j . Using it we can 

define the Gaufi curvature K q = det(<7 s ? ) and the 
mean curvature H q = Tr (g s q ). 

In the special case of plane curves (M = S 1 and 
d = 2) we use the letter c for the immersed curve and 
we denote the corresponding unit length tangent vector 
by 

C , 2 

v = v c = — r = —Jn c where J = y—1 on C = K . 

The covariant derivative of the pull-back metric reduces 
to arclength derivative, and the induced volume form 
to arclength integration: 



D„ = 



ds = \c e \d0. 



Using this notation the length of a curve can be written 
as 



4 = 



ds . 



s 1 



In this case GauB and mean curvature are the same and 
are denoted by k = (D s v,n). 

2.2 Riemannian submersions 

In this article we will repeatedly induce a Riemannian 
metric on a quotient space using a given metric on the 
top space. The concept of a Riemannian submersion will 
allow us to achieve this goal in an elegant manner. We 
will now explain in general terms what a Riemannian 
submersion is and how geodesies in the quotient space 
correspond to horizontal geodesies in the top space. 

Let (E, Ge) be a possibly infinite dimensional weak 
Riemannian manifold; weak means that Ge '■ TE — > 
T* E is injective, but need not be surjective. Let Q be 
a smooth possibly infinite dimensional Lie group and 
Q x E — > E be a smooth right group action on E and 
assume that B := E/Q is a manifold. Let 7r : E — > B be 
the projection, which then is a submersion of smooth 
manifolds, that is, Ttt : TE —> TB is surjective. Then 

V = V(tt) := ker(TTr) C TE 

is called the vertical subbundle. Assume that Ge is in 
addition invariant under the action of Q . Then the pro- 
jection 7r : E — > B is a Riemannian submersion: for 
x G E and Y € T^i x \B the expression 



defines a weak Riemannian metric Gb on B. Some- 
times the the G^-orthogonal space V^) C TE is a 
fiber linear complement in TE. In general, the orthog- 
onal space is a complement (for the closure of V(tt)) 
only if taken in the fiberwise G^-completion TE of TE; 
see |84| for the notion of a robust Riemannian manifold 
which takes care of the existence of this fiberwise com- 
pletion as a smooth bundle over E, and of the existence 
of the geodesic equation. In section[8]we have to use TE. 
The horizontal subbundle Hor = Hor(7r, G) is the Ge~ 
orthogonal complement in TE or in TE, respectively. 
Then any vector X € TE can be decomposed uniquely 
in vertical and horizontal components as 



X = X" 



X 



hor 



and the mapping 

T x Tr\nor x ■ Hor a -> T n ^B or T^B 

is a linear isometry of (pre) Hilbert spaces for all x € E. 

Theorem 2.1 Consider a Riemannian submersion n : 
E — > B, and let 7 : [0, 1] — > E be a geodesic in E. 

1. If 7' (t) is horizontal at one t, then it is horizontal 
at all t. 

2. If j'(t) is horizontal then tt o 7 is a geodesic in B. 

3. If every curve in B can be lifted to a horizontal curve 
in E, then there is a one-to-one correspondence be- 
tween curves in B and horizontal curves in E. This 
implies that instead of solving the geodesic equation 
on B one can equivalently solve the equation for hor- 
izontal geodesies in E. 



See 87 Sect. 26] for a proof, and 84 for the case of 



robust Riemannian manifolds. 

3 The spaces of interest 

3.1 Immersions and embeddings 

Parametrized surfaces will be modeled as immersions 
or embeddings of the configuration manifold M into 
M. d . We call immersions and embeddings parametrized 
since a change in their parametrization (i.e., applying a 
diffcomorphism on the domain of the function) results 
in a different object. We will deal with the following 
sets of functions: 



Emb(M, R d ) C Imnij (M, R d ) 



C Imm(M,K d ) C C°°(M,K d ) . 



(2) 



\y\\l E 



inf{||X| 



X e T X E, Tn.X = Y) 



Here C°°(M,R d ) is the set of smooth functions from M 
to R d , Imm(M,]R ci ) is the set of all immersions of M 
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into R d , i.e., all functions q € C°°(M, E d ) such that T^g 
is injective for all x € M. The set Imm/(M, M d ) consists 
of all free immersions q; i.e., the diffcomorphism group 
of M acts freely on q, i.e., qo ip = f implies ip = Id^ for 
all ip 6 Diff(M). By (25j Lem. 3.1], the isotropy group 
Diff(M), := € Diff(M) : qo/ = g} of a general 
immersion q is always a finite group which acts strictly 
discontinuously on M so that M — > M/T>iS(M) q is a 
covering map. Emb(M , N) is the set of all embeddings 
of M into R d , i.e., all immersions q that are homeomor- 
phisms onto their image. 

Theorem 3.1 The spaces Imm(M, R d ), Imm/(M, R d ) 
and Emb(Af, R d ) are Frechet manifolds. 

Proof Since M is compact by assumption it follows that 



C°°(M,R d ) is a Frechet manifold 72 Sec. 42.3]. All 



inclusions in ^ are inclusions of open subsets: first 
Imm(M, R d ) is open in C°°(M,R d ) since the condition 
that the differential is injective at every point is an open 
condition on the one- jet of q |85| Sec. 5.1]. Imm/(M,IR d ) 
is open in Imm(M,K d ) by [25[_Thm. 1.5]. Emb(Af,M d ) 
is open in lmm f (M,R d ) by [72] Thru. 44.1]. Thus all 
the spaces are Frechet manifolds as well. 



3.2 Shape space 

Unparametrized surfaces are equivalence classes of pa- 
rametrized surfaces under the action of the reparame- 
trization group. 



Theorem 3.2 (Thm. 1.5, |25| ) The quotient space 
B hf (M,R d ) := Imm/(M,M d )/DifF(M) 
is a smooth Hausdorff manifold and the projection 
n : rmm/(M,R d ) -> Bij{M,R d ) 

is a smooth principal fibration with Diff (M) as structure 
group. 

Forqe lmm f (M,R d ) we can define a chart around 
tt(?) G Bij(M, R d ) by 

noil; q : C°°(M, (-e,e)) -)■ B tJ (M,R d ) 

with e sufficiently small, where 

ip q : C°°(M, (-£,£)) -> lmm f (M,R d ) 

is defined by ip q (a) = q + an q and n q is the unit-length 
normal vector to q. 

Corollary 3.3 The statement of Thm. \3.S\ does not 
change, if we replace Bij(M,R d ) by B e (M,R d ) and 
Imm(M,R d ) byEmb(M,R d ). 



As Emb(Af, R d ) is an open subset of Imm / (M,R d ) 
and is Diff (M )-equi variant, the quotient 

B e (M, R d ) := Emb(A/,]R' i )/Diff(A//) 

is an open subset of Bij(M,R d ) and as such itself a 
smooth principal bundle with structure group Diff (M). 

3.3 Some words on orbifolds 
The projection 

Imm(Af,R d ) -> Bi(M,R d ) := Imm(Af, R d )/ Diff(M) 

is the prototype of a Riemannian submersion onto an 
infinite dimensional Riemannian orbifold. In the arti- 
cle |115[ Prop. 2.1] it is stated that the finite dimen- 
sional Riemannian orbifolds are exactly of the form 
Al/G for a Riemannian manifold M and a compact 
group G of isometries with finite isotropy groups. Cur- 
vature on Riemannian orbifolds is well defined, and it 
suffices to treat it on the dense regular subset. In our 
case Bij(M, R d ) is the regular stratum of the orbifold 
Bi(M, R d ). For the behavior of geodesies on Rieman- 
nian orbit spaces M/G see for example [l]; the easiest 
way to carry these results over to infinite dimensions 
is by using Gauss' lemma, which only holds if the Rie- 
mannian exponential mapping is a diffeomorphism on 
an Gi mm -open neighborhood of in each tangent space. 
This is rarely true. 

Given a Diff (Af )-invariant Riemannian metric on 
Imm(M, R d ), one can define a metric distance dist Si on 
Bi(M,R d ) by taking as distance between two shapes 
the infimum of the lengths of all (equivalently, hor- 
izontal) smooth curves connecting the corresponding 
Diff (AZ)-orbits. There are the following questions: 

— Does dist Bl separate points? In many cases this has 
been decided. 

— Is (Bi(M, K d ),dist Si ) a geodesic metric space? In 
other words, does there exists a rectifiable curve 
connecting two shapes in the same connected com- 
ponent whose length is exactly the distance? This is 
widely open, but it is settled as soon as local mini- 
mality of geodesies in Imm(Af, R d ) is established. 

3.4 Diffeomorphism group 

Concerning the Lie group structure of the diffeomor- 
phism group we have the following theorem. 

Theorem 3.4 (Thm. 43.1, (73)) Let M be a smooth 
manifold, not necessarily compact. The group 

Diff c (M) = {ip : ip,ip- x G C°°(M,M), 

supp(</? — Id) compact} 
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of all compactly supported diffeomorphisms is an open 
submanifold of C°° [M, M) and composition and inver- 
sion are smooth maps. It is a regular Lie group and the 
Lie algebra is the space X C (M) of all compactly sup- 
ported vector fields. 

An infinite dimensional smooth Lie group G with 
Lie algebra g is called regular, if the following two con- 
ditions hold: 

— For each smooth curve X £ C ca (M.,g) there exists 
a unique smooth curve g £ C°°(R,G) whose right 
logarithmic derivative is A, i.e., 



5(0) = e 

d t g(t) =WW)X(t) = X(t).g(t) 



(3) 



— The map evolg(X) = g(l) where g is the unique 
solution of (J3j , considered as a map between the 
spaces evoiG : C°°(]R,fj) -> G is C°° -smooth. 

If M is compact, then all diffeomorphisms have com- 
pact support and Diff c (M) = Diff(M). For R n the 
group Diff(]R d ) of all orientation preserving diffeomor- 
phisms is not an open subset of C°°(R , R d ) endowed 
with the compact C°°-topology and thus it is not a 
smooth manifold with charts in the usual sense. There- 
fore, it is necessary to work with the smaller space 
Diff c (M d ) of compactly supported diffeomorphisms. In 
Sect. u\ we will also introduce the groups Difff/oo (R d ) 
and Diff^R' 1 ) with weaker decay conditions towards 
infinity. Like Diff c (R d ) they are smooth regular Lie 
groups. 

3.5 The space of Riemannian metrics 

We denote by Met(M) the space of all smooth Rieman- 
nian metrics. Each g £ Met(M) is a symmetric, positive 
definite (°) tensor field on M, or equivalently a point- 
wise positive definite section of the bundle S 2 T*M. 



Theorem 3.5 (Sec. 1.1, |50|) Let M be a compact 
manifold without boundary. The space Met(M) of all 
Riemannian metrics on M is an open subset of the 
space r \S 1 T* M) of all symmetric (°) tensor fields and 
thus itself a smooth Frechet-manifold. 

For each g £ Met(M) and x £ M we can regard 
g(x) as either a map 

g{x) : T X M x T X M ->■ R 

or as an invertible map 

g x :T x M^T*M. 

The latter interpretation allows us to compose g, h £ 
Met(M) to obtain a fiber-linear map g~ 1 -h : TM — > 
TM. 



4 The i 2 -metric on plane curves 

4.1 Properties of the L 2 -metric 

We first look at the simplest shape space, the space of 
plane curves. In order to induce a metric on the mani- 
fold of un-parametrized curves Bij(S 1 , M 2 ) we need to 
define a metric on parametrized curves Imm(5 1 , R 2 ), 
that is invariant under reparametrizations, c.f. Sect. |2.2j 
The simplest such metric on the space of immersed 
plane curves is the L 2 -type metric 

G° c (h,k) = { (h(9),k(6))ds. 
Js 1 

For this metric the horizontal bundle consists of all tan- 
gent vectors h that are pointwise orthogonal to c', i.e., 
h(9) = a(0)n c {9) for some scalar function a £ C ca (S 1 ). 
Using the charts from Thm. |3.2| the metric on the quo- 
tient space B i j(S 1 ,R 2 ) thus reduces to 

G° c (a,b) = [ a(6)b(d)ds. 



This metric was first studied in the context of shape 



analysis in 89 . The geodesic equation for the G°-metric 
on Imm/(5' 1 ,R 2 ) is given by 



By Thm. 2.1 geodesies on B i j(S 1 ,M 2 ) correspond to 
horizontal geodesies on Imm(S' 1 ,R 2 ), which are of the 
form ct = a.n Cl with a scalar function a(i, 9). Thus the 
geodesic equation Q reduces to an equation for a(t, 9), 



a, 



-na 



Note that this is not an ODE for a, because k c , being 
the curvature of c, depends implicitly on a. It is however 
possible to eliminate n and arrive at 



a-tt 



a 6 a09 a 6 a 9 e 
2v 4 v 5 



a(0,9) 







v{9) 



VMo,0)| ' 

a nonlinear hyperbolic PDE of second order. 

Open question. Are the geodesic equations on either 
of the spaces Imm(S' 1 ,]R 2 ) or B i j(S 1 ,R 2 ) for the L 2 - 
metric (locally) well-posed? 

The L 2 -metric is among the few for which the sec- 
tional curvature on B i j(S 1 ,R 2 ) has a simple expres- 
sion. Let C = 7r(c) £ B^fi^S 1 ,R 2 ) and choose c £ 
Imm/(S' 1 ,R 2 ) such that it is parametrized by constant 
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speed. Take a.n c , b.n c G Horgo(c) two orthonormal hor- 
izontal tangent vectors at c. Then the sectional curva- 
ture of the plane spanned by a, b G T c B i j(S 1 ,M. 2 ) is 
given by the Wronskian 



kc{P(a,b)) 



(abg — agb) ds . 



(5) 



In particular the sectional curvature is non-negative 
and unbounded. 

Remark 4-1 This metric has a natural generalization 
to the space Imm(M, M. d ) of immersions of an arbitrary 
compact manifold M. This can be done by replacing 
the integration over arc-length with integration over the 
volume form of the induced pull-back metric. For q G 
Imm(M, M d ) the metric is defined by 

G°(M)= / (h(x),k(x))vol(q*g). 
Jm 



The geodesic spray of this metric was computed in 17 
and the curvature in [62] . 

For all its simplicity the main drawback of the L 2 - 
metric is that the induced geodesic distance vanishes. 
If c : [0, 1] -> Imm(S' 1 ,R 2 ) is a path, denote by 



Len T 



Imm(c) = / yjG c(t) (ct(t),Ct(t))dt 



its length. The geodesic distance between two points is 
defined as the infimum of the pathlength over all paths 
connecting the two points, 

dist imm (co , ci ) = inf Lenf mm (c) . 

c(0)=c o , 
c(l)=ci 

For a finite dimensional Riemannian manifold (M, G) 
this distance is always positive, due to the local invert- 
ibility of the exponential map. This does not need to 
be true for a weak Riemannian metric in infinite dimen- 
sions and the i 2 -metric was the first known example, 
where this was indeed false. We have the following re- 
sult. 

Theorem 4.2 The geodesic distance function induced 
by the metric G° vanishes identically on Imm(5' 1 ,K 2 ) 
and Bij{S\R 2 ). 

For any two curves Co, ci G Imm(5 1 ,K 2 ) and e > 
there exists a smooth path c : [0, 1] —> Imm(S' 1 , M 2 ) with 
c(0) = Co, c(l) = c\ and length Len^ nm (c) < e. 

For the space B i j(S 1 ,M. 2 ) an explicit construction 



of the path with short length was given in 89 . Heuris 



tically, if the curve is made to zig-zag wildly, then the 
normal component of the motion will be inversely pro- 
portional to the length of the curve. Since the normal 



component is squared the length of the path can be 
made arbitrary small. This construction is visualized in 
FigH 

For Imm( 1 5' 1 , M 2 ) vanishing geodesic distance is pro- 
ven in [7], and it is a combination of the results on 
BijiS 1 ^ 2 ) and on Diff(5' 1 ). 




0.2 0.3 



0.6 0.7 0.8 0.9 



Fig. 2 A short path in the space of un-parametrized curves 
that connects the bottom to the top line. Original image pub- 
lished in [89]. 



Remark 4-3 In fact, this result holds more generally for 
the space Imm(M,K d ). One can also replace M. d by an 
arbitrary Riemannian manifold N; see 



The vanishing of the geodesic distance leads us to 
consider stronger metrics that prevent this behavior. 
In this article we will present three different classes of 
metrics: 

— Almost local metrics: 

G*(h,k)= [ *(q)(h,k)vol(q*g), 



where ^ : Imm(M,M d ) 
function. 

Sobolev type metrics: 



is a suitable smooth 



G L q {h,k)= / (L q h,k)vol(q*g), 
Jm 

where L q : T g lmm(M,R d ) -> 
suitable differential operator. 

Metrics that are induced by right invariant metrics 
on the diffeomorphism group of the ambient space. 



T q lmm(M,R d ) is 
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4.2 Gradient flows on curves 

The L 2 -metric is used in geometric active contour mod- 
els to define gradient flows for various energies. For ex- 
ample the curve shortening flow 

is the gradient flow of the energy E(c) — j sl ds = l c 
with respect to the L 2 -metric. 

The following example is taken from [81] . The cen- 
troid based energy E(c) = ||/u(c) — w\ 2 , with w £ R 2 
fixed and /i(c) = j- J s± cds denoting the centroid, at- 
tains its minimum when ^t(c) = w. The L 2 -gradient of 
the energy is 

V°-E(c) = (fi(c) - w, n c )n c + K c (/i(c) - c, (j,(c) - w)n c . 
We see from the second term that the gradient flow 
ct = -V°E(c) 

tries to decrease the length of the curve for points with 
(n(c) — c, n{c) — w) < 0, but increase for (11(c) — c, /i(c) — 
w) > 0. This latter part is ill-posed. However the ill- 
posedness of the gradient flow is not an intrinsic prop- 
erty of the energy, it is a consequence of the metric we 
chose to define the gradient. For example the gradient 
flow with respect to the 7J 1 -metric 

Gl(h,k)= [ ±(h,k)+£ c (D s h,D s k)ds 
Js 1 ' 



is locally well-posed. See 116 118| for more details on 
Sobolev active contours and applications to segmenta- 
tion and tracking. The same idea has been employed 



for gradient flows of surfaces in 131 



5 Almost local metrics on shape space 



Almost local metrics are metrics of the form 



M 



$(Vol q , H q ,K q )(h, k)vol(q*g) 



where <P <E G°°(]R 3 , M >0 ) is a function of the total vol- 
ume Vol g , the mean curvature H q and the GauB cur- 
vature K q . The name "almost local" is derived from 
the fact that while H q and K q are local quantities, the 
total volume Vol g induces a mild non-locality in the 
metric. If <P — <£(Vol) depends only on the total vol- 
ume, the resulting metric is conformally equivalent to 
the L 2 -metric, the latter corresponding to <P = 1. 
More generally almost local metrics are defined as 

0*(h,k)= f &(q)(h,k)vol(q*g), 



where <f : Imm(M,R d ) -> C o °(M,R >0 ) is a smooth 
function that is equivariant with respect to the action 
of Diff(M), i.e., 

<P(q oip) = &(q) oip, qe Imm(Af , R d ) , ip G Diff(M) . 

Equivariance of ^ then implies the invariance of G* and 
thus G* induces a Riemannian metric on the quotient 
B tJ (M,R d ). 

For an almost local metric G* the horizontal bundle 
at q S Imm/(Af , R d ) consists of those tangent vectors 
h that are pointwise orthogonal to q, 

Hor*^) = {he T q Imm f (M,M. d ) : (h(x),n q (x)) =0} . 



Using the charts from Thm. 3.2 the metric G* on 
B zJ (M,R d ) is given by 



Gt (q) (a,b)= f V(q)abvo\(q*g) : 

J M 



with a, be T^^B.jiM.R" 1 ). 

Almost local metrics, that were studied in more de- 
tail include the curvature weighted G A -metrics 



Gf(h,k) = / (1 + An 2 c )(h,k)ds. 
Js 1 



(6) 



with A > in 89 and the conformal rescalings of the 
L 2 -metric 

Gt(h,k) = $(Q f (h,k)ds, 



with <Z> e G oo (M >0 , K> ) in [l09][T26] , both on the space 
of plane curves. More general almost local metrics on 
the space of plane curves were considered in [90] and 
they have been generalized to hypersurfaces in higher 
dimensions in [3 13 14 . 



5.1 Geodesic distance 

Under certain conditions on the function & almost local 
metrics are strong enough to induce a point-separating 
geodesic distance function on the shape space. 

Theorem 5.1 If & satisfies one of the following con- 
ditions 

1. ■%) > l + AH 2 

2. <%) > A\o\ q 

for some A > 0, then the metric G* induces a point- 
separating geodesic distance function on Bij(M,R d ) , 
i.e., for Go 7^ C\ we have distfj. ^(GojGi) > 0. 
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For planar curves the result under assumption [T] is 
proven in |89[ Sec. 3.4] and under assumption [2] in |109| 
Thm. 3.1]. The proof was generalized to the space of 
hypersurfaces in higher dimensions in [l3j Thm. 8.7]. 

The proof is based on the observation that under 
the above assumptions the G^-length of a path of im- 
mersions can be bounded from below by the area swept 
out by the path. A second ingredient in the proof is the 
Lipschitz-continuity of the function ^/Vol 9 . 

Theorem 5.2 If W satisfies 

${q) > 1 + AH 2 , 

then the geodesic distance satisfies 



x/Vol^- v/VoI^ 



< 



2V A 



distg (Qi,Q 2 



In particular the map 

VWl : (B il/ (Af,R <, ),disti* / ) -> M> 

is Lipschitz continuous. 

This result is proven in [89j Sec. 3.3] for plane curves 
and in [13[ Lem. 8.4] for hypersurfaces in higher dimen- 
sions. 

In the case of planar curves [109] showed that for 
the almost local metric with &(c) — l c the geodesic 
distance on Bij(S 1 ,M. 2 ) is not only bounded by but 
equal to the infimum over the area swept out, 



distil (Co, CO 



inf / |det dc(i,0) | d0df. 

7r(c(0))=C ./six fo.i] 



Remark 5.3 No almost local metric can induce a point 
separating geodesic distance function on Immf (M, R d ). 
When we restrict the metric G* to an orbit qo Diff(M) 
of the Diff (M)-action, the induced metric on the space 
goDiff(M) = Diff(M) is a right-invariant weighted L 2 - 
type metric, for which the geodesic distance vanishes. 
Thus 

^stf mm (q,qoif) = 0, 

for all q G Imm/(M,R d ) and ip G Diff(M). See Sect. 
17.41 for further details. 



It is possible to compare the geodesic distance on 
shape space with the Frechet distance. The Frechet dis- 
tance is defined as 



dists, ,{Qo,Qi) 



inf || g 

go, si 



(7) 



where the infimum is taken over all immersions <?o,9i 
with n(qi) — Qi. Depending on the behavior of the 
metric under scaling, it may or may not be possible to 
bound the Frechet distance by the geodesic distance. 



Theorem 5.4 (Thm. 8.9, [l3]) IfW satisfies one of 
the conditions, 

1. *{q) < Ci + C 2 Hf 

2. *{q)<C x Vo\ k q , 

with some constants C\,C-2 > and k < , then 
there exists no constant C > 0, such that 

dist|~ (Qo, Qi)<C dist* ^ (Qo, Qi) , 

holds for all Q ,Qi G B i>f (M,M. d ). 

Note that this theorem also applies to the G^-metric 
for planar curves defined in ([6]). Even though Thm. 5.4 
states that the identity map 



),dist 



G 



, dist 1 



is not Lipschitz continuous, it can be shown that is 
continuous and thus the topology induces by dist is 
stronger than that induced by dist L . 

Theorem 5.5 (Cor. 3.6, |89|) The identity map on 
Bij(S l ,M. 2 ) is continuous from {B. lJ (S 1 ,R 2 ),dist GA ) 
to {Bij{S 1 , K 2 ), dist L ) and uniformly continuous on 
every subset, where the length lc is bounded. 

As a corollary to this result we obtain another proof 
that the geodesic distance for the G^-metric is point- 
separating on Bij(S 1 ,M. 2 ). 



5.2 Geodesic equation 

Since geodesies on Bij correspond to horizontal geo- 
desies on Imm/(M, K d ), see Thm. 2.1 to compute the 
geodesic equation on Bij(M, M. d ) it is enough to re- 
strict the geodesic equation on Immf(M,Mr) to hori- 
zontal curves. 

As an example for the resulting equations we will 
present the geodesic equations on Bij(M, R d ) for the 
almost local metric with ^(q) = 1 + AH 2 , which is 
a generalization of the metric ([6]), and the family of 
metric &(q) = $(Vol g ), which are conformal rescalings 
of the L 2 -metric. 

Theorem 5.6 Geodesies of the almost local G 1, -metric 
with &(q) = 1 + AH 2 on B it f(M,R d ) are given by so- 
lutions of 



qt = an q , 

1„2 rr2 



9 = q*g , 



2 - a H q 



zgfejgiHgWa + 2aV q H q , V 9 a) 
^^H q + Tr((g-h q ) 2 )) . 



A. 
l+AHl 
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For the family of metrics with ^(q) 
are given by 

b(t) 



<2>(Volg) geodesies 



It = 



at 



<2>(Vol Q 



g = q 9, 

<P'(Vol q ) 



a 2 vol(ff) 



For the G^-metric and planar curves the geodesic equa- 



tion was calculated in 89 Sect. 4.1], whereas for con- 



formal metrics on planar curves it is presented in 109 



Sect. 4]. For hypersurfaces in higher dimensions the 



equations are calculated in 13 Sect. 10.2 and 10.3]. 

Note that both for A = and 'P(q) = 1 one recovers 
the geodesic equation for the L 2 -metric, 

a t = \H 2 q a 2 . 

Similarly to the case of the L 2 -metric it is unknown, 
whether the geodesic equations are well-posed. 

Open question. Are the geodesic equations on either 
of the spaces Imm(M, R 2 ) or Bij(M, E 2 ) for the almost 
local metrics (locally) well-posed? 



is constant along geodesies. 

If the function W satisfies the scaling property 



&(\q) = X 



_ \-dim(M)-2 



W{q), q E Imm(Af,R d ), A G K >0 , 



then the induced metric G* is scale invariant. In this 
case the scaling momenta are conserved along geodesies 
as well: 



&(q)(q,qt)voi(g) 



scaling momentum 



M 



For plane curves the momenta are 

/ &(c)(c' ,Ct)fj,ds reparametrization momentum 
Js 1 



S 1 



Si 



<P(c)ct ds 
(Jc,c t ) ds 



&(c)(c,ct)di 



linear momentum 
angular momentum 
scaling momentum 



with /i £ X(S 1 ) and Jc denoting rotation by ?. 



5.3 Conserved quantities 

If the map & is equivariant with respect to the DifF (Ab- 
action, then the G*-metric is invariant, and we obtain 
by Noether's theorem that the reparametrization mo- 
mentum is constant along each geodesic. The repara- 
metrization for the G^-metric is given by 

^(q)9(qj-) ™\(q*g) € F(T*M ® M vol(M)) , 

with g = q*g and the pointwise decomposition of the 
tangent vector q t — qj + qf- of q t into q^~ = g(qt,n q )n q 
and qj (x) € T q t x \q(M). This means that for each X G 
3C(M) we have 

/ &(q)g(qJ,X) vo\(q*g) = const. 

J M 

If & is additionally invariant under the action of 
the Euclidean motion group R d x SO(d), i.e., &(0.q + 
v) = \P(q), then so is the G*-metric and by Noether's 
theorem the linear and angular momenta are constant 
along geodesies. These are given by 

f $(q)q t vol(q*g) € R d 

J M 

[ V(q)qAq t vol(q*g) G f\ R d ^ so(d)* . 
J M 

The latter means that for each fi G so(d) the quantity 

/ &{q)g{&-q> qt) vo\(q*g) = const. 

J M 



5.4 Completeness 

Regarding geodesic completeness, one can look at the 
set of spheres with a common center. This set is a to- 
tally geodesic submanifold of Bij{S 1 , M. d ). One can 
explicitly calculate the length of a geodesic as spheres 
shrink towards a point and when they expand towards 
infinity. When it is possible to shrink to a point with 
a geodesic of finite length, the space can obviously not 
be geodesically complete. This is the case under the 
following conditions. 

Theorem 5.7 (Thm. 9.1, [l3]) If\P satisfies one of 
the conditions, 

1. $(q) < Gi + C 2 H 2k 

2. $(z)<CxY6§, 

with some constants C\,C2 > and k < — i*^, then 
the spaces lmm(S d ~ 1 , R d ) and B iJ (S d - 1 , R d ) are not 
geodesically complete with respect to the G* -metric. 

Note that these arc the same conditions as in Thm. 
|5.4| For other choices of M scalings will in general not 
be geodesic, but under the same condition an immersion 
can be scaled down to a point with finite energy. What 
conditions are sufficient to prevent geodesies from de- 
veloping singularities and thus make the spaces geodesi- 
cally complete is unknown. 

Concerning metric completeness, it cannot be ex- 
pected that a weighted L 2 -type metric will be able to 
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prevent immersions from losing smoothness in the com- 
pletion. We have only a partial result available for the 
G^-metric ^ on plane curves. 

Similarly to the definition of B it f(S l , K 2 ), we can 
define the larger space 

£f p (S\R 2 ) = Lip(S\R 2 )/~ 

of equivalence classes of Lipschitz curves. We identify 
two Lipschitz curves, if they differ by a monotone cor- 
respondence. This can be thought of as a generalization 
of reparametrizations, which allow for jumps and inter- 
vals of zero speed; see [89[ Sect. 2.11]. Equipped with 
the Frechet-distance Q, the space Si tp (S rl ,M 2 ) is met- 
rically complete. 

A curve G G B? p (S 1 ,R 2 ) is called a 1-BV rectifiable 
curve, if the turning angle function a of an arc-length 
parametrized lift c £ Lip(S' ,M ) of G is a function of 
bounded variation. 

Theorem 5.8 (Thm. 3.11, (89]) The completion of 
the metric space (Bij(S 1 ,M 2 ),dist G ) is contained in 
the shape space i?| lp (5 1 ,K 2 ) of Lipschitz curves and it 
contains all 1-BV rectifiable curves. 



5.5 Curvature 

The main challenge in computing the curvature for al- 
most local metrics on Imm(M, M. d ) is finding enough 
paper to finish the calculations. It is probably due to 
this that apart from the L 2 -metric we are not aware of 
any curvature calculations on the space Imm(M, R d ). 
For the quotient space Bij(M, M. d ) the situation is a bit 
better and the formulas a bit shorter. For plane curves 
and conformal metrics the curvature has been calcu- 
lated in [109] and for t^(c) = $(£ c ,k c ) in f9Cfl. Similar 



for higher dimensional surfaces the curvature has been 
calculated for \P(q) = <P(Vol q , H q ) in (13 . 



The sectional curvature for the L 2 -metric on plane 
curves ^ is non-negative. In general the expression for 
the sectional curvature for almost local metrics with 
<P ^ 1 will contain both positive definite, negative def- 
inite and indefinite terms. For example the sectional 



curvature of the metric (|6j). with <^(c) 
plane curves has the following form. 



= 1 + An 2 



Theorem 5.9 (Sec. 4.6, [89]) Let C £ B i j(S 1 ,R 2 ) 
and choose c £ Imrn.(S' 1 ,R 2 ) such that C = 7f(c) and c 
is parametrized by constant speed. Let a.n, b.n £ Hor(c) 
two orthonormal horizontal tangent vectors at c. Then 
the sectional curvature of the plane spanned by a,b £ 



T c B iJ (S 1 ,R 2 ) for the G A metric is 

k c (P(a, &)) = -/ A{a.D 2 s b - D 2 s a.b) 2 ds 
Js 1 

(1 - An 2 ) 2 - AA 2 k.D 2 k + 8A 2 (D s k) 



2(1 + Ak 2 ) 



(a.D s b- D s a.b) 2 ds 



It is assumed, although not proven at the moment, 
that for a generic immersion, similar to Thm. |7.14[ the 
sectional curvature will assume both signs. 



5.6 Examples 

To conclude the section we want to present some ex- 
amples of discrete solutions to the geodesic boundary 
value problem for given shapes Qo,Qi £ Bij(M,M. d ). 
One method to tackle this problem is to directly mini- 
mize the horizontal path energy 



M 



${Vo\ q ,H q )(q u n q ) 2 vo\{q*g) 



over the set of paths q of immersions with fixed end- 
points qo,qi that project onto the target surfaces C/o 
and Qi, i.e., n(qi) — Qi. The main advantage of this 
approach for the class of almost local metrics lies in 
the simple form of the horizontal bundle. Although we 
will only show one specific example in this article it is 
worth to note that several numerical experiments are 
available; see: 



90 



126 



for the G^-metric and planar curves. 



127) for conformal metrics and planar curves. 
3 13 ,14 for surfaces in M 3 . 




Fig. 3 A geodesic in the G^-metric joining two shapes of 
size about 1 at distance 5 apart with A = .25, using 20 time 
samples and a 48-gon approximation for all curves. Original 
image published in |89|. 



The example we want to present here, is concerned 
with the behavior of the G^-metric matching curves 



that are far apart in space. In the article 89 the au- 



thors showed that pure translation of a cigar-like shape 
with a cross-section of 2\[A. is (locally) a geodesic for 
the G A -metric. Thus one might expect that a geodesic 
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between distant curves will asymptotically utilize this 
cigar shaped curve, translate this optimal curve and 
then deform it to the target shape. In fact the numerical 
examples resemble this behavior as can be seen in Figj3j 
Note that the cross-section of the middle figure - which 
is highlighted - is slightly bigger than 2yfA. A reason 
for this might be that the distance between the two 
boundary shapes is not sufficiently large. In the arti- 
cle [13] it has been shown that this behavior carries over 
to the case of higher-dimensional surfaces, c.f. Fig. [4j 
Note that the behavior of the geodesies changes dra- 
matically if one increases the distance further, namely 
for shapes that are sufficiently far apart the geodesies 
will go through a shrink and grow behavior. This phe- 
nomenon is based on the fact that it is possible to shrink 
a sphere to zero in finite time for the G^-metric. Then 
geodesies of very long translations will go via a strong 
shrinking part and growing part, and almost all of the 
translation will be done with the shrunken version of 
the shape. This behavior, which also occurs for the class 
of conformal metrics, is described in 
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Fig. 4 Middle figure of a geodesic between two unit spheres 
at distance 3 apart for A = 0.2, A = 0.4, A = 0.6, A = 0.8. In 
each of the simulations 20 time steps and a triangulation with 
720 triangles were used. Original image published in [131. 



6 Sobolev type metrics on shape space 

Sobolev-type inner metrics on the space Imm(Af, M. d ) 
of immersions are metrics of the form 



G L Ah,k) 



(L q h,k)vol(q*g) 



M 



where for each q £ M , L q is a pseudo-differential oper- 
ator on T g Imm(M, M. d ). To be precise we assume that 

L : TImm(M,M d ) ->■ TImm(M,R d ) 

is a smooth base-point preserving bundle isomorphism, 
such that for every q 6 Imm(M, M. d ) the map 



L q : T q Imm(Af, ] 



T q Imm(Af, R d ) 



is a pseudo-differential operator, that is symmetric and 
positive with respect to the L 2 -metric. An example for 
such an operator L is 



where A 9 is the Laplacian of the induced metric g = q* g 
on M. 

We will also assume that L is invariant under the 
action of the reparametrization group Diff(A/), i.e., 

(Lqh) o ip — Lq v (ho If) , (9) 

for all <f, q and h. Then the metric G L is invariant under 
Diff(Af) and it induces a Riemannian metric on the 
quotient space Bij(M,R d ). 

In contrast to the class of almost local metrics, for 
whom the horizontal bundle of the submersion 

Imm/(M,R d ) -> Bij(M,R d ) 

consisted of tangent vectors, that are pointwise orthog- 
onal to the surface, here the horizontal bundle cannot 
be described explicitly. Instead we have 

Ror L (q) = T q lmm f (M,R d ) : L q = an q } , 

where a € C°°(M, R) is a smooth function. Thus to 
parametrize the horizontal bundle we need to invert 
the operator L q . 

General Sobolev-type inner metrics on the space of 
immersed plane curves have been studied in 



90 and 



on surfaces in higher dimensions in |11| . Numerical ex- 
periments for special cases of order one Sobolev type 



metrics are presented in the articles 4 59 114 




0. 180 



0.154 



0.129 



0.103 



0.O77J 



0.O5 14 



0.0257 



0.00 



L q h = h+(A 9 ) l h, l>0, 



(8) 



Fig. 5 In this figure we show the Karcher mean of five vase- 
shaped objects with respect to the Sobolev metric of order 
one - as defined in (|8b - on the space of parametrized surfaces 
Imm(5 1 x [0, 1], R 3 )TThe mean shape, which is displayed in 
the center of the figure is computed using an iterated shooting 
method. The colored regions on the averaged shapes encode 
the Euclidean length of the initial velocity of the geodesic, 
which connects each shape to the mean. The color of the 
mean was chosen for artistic purposes only. Original image 
published in [41. 



In [15] the authors consider metrics of the form 
G L q {h, k) = #(Vol 9 ) / (L q h, k) vo\(q*g) . 
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These are a combination of Sobolev-type metrics with 
a non-local weight function, that can be chosen such 
that the resulting metric is scale-invariant. Sobolev type 
metrics are far better investigated and understood on 
the the manifold of plane curves then in higher dimen- 
sion. Therefore, we will discuss this case separately in 

sect, urn 



The image of the space ImmfS 1 , M 2 )/ transl. under 
the R-transform is a co-dimension 2 submanifold of the 
fiat space C°°(S 1 ,M. 2 ). 

This representation not only allows to efficiently dis- 
cretize the geodesic equation, but also to compute the 
curvature of Imm(S ll ,R 2 ); See 5 114 for details. 



6.1 Sobolev metrics on plane curves 

A reparametrization invariant Sobolev-type metric on 
the space of pl ane curves has been first introduced by 
Younes in 
metric 



128 



There he studied the homogeneous H 1 



G?\h,k) = ( (D s h,D s k)ds. 
Js 1 



However, this is not a metric on Imm(S' 1 ,IR 2 ) but only 
on the quotient space Irnm(S' 1 ,R 2 )/ transl. In order to 
penalize bending and stretching of the curve differently 
it has been generalized in 94 95 to 

G a /(h,k)- 



a-(D s h,n c )(D s k,n c ) 

+ b 2 (D s h 1 v c )(D s k,v c ) ds 



In this metric the parameters a, b can be interpreted as 
the tension and rigidity coefficients of the curves. For 
a = 1,6 = i a computationally efficient representation 
of this metric - called the Square Root Velocity Trans- 
form (SRVT) - has been found in |114| and it has been 
generalized for arbitrary parameters a, b in [5l. Follow- 



114 we will describe this transformation for the 



mg 

case a = \ , b = 1 : 



R . |lmm(S' 1 , 



R 2 )/ transl. -> C°°(S' 1 ,K 2 ) 
c H> y^fc/fw. 

The inverse of this map is given by 



R- 



C°°(S' 1 ,]R 2 ) -> Imm([0,27r],IR 2 )/ transl. 
e ^ c + J°\e(8M6)d9. 



Note that R 1 (e) is a closed curve if and only if 



i{0)\e{0) d0 = O. 



Theorem 6.1 Consider the fiat 1? -metric 



G 



ZAflat 



(e,f) 



s 1 



(ej) d6 



on C°°{S\R 2 ). The pullback of the metric G i2 > flat by 
the R-transform is the elastic metric with coefficients 
a = 1, b = jj. 









Fig. 6 A geodesic in the shape space £?, f(S 1 ,M. 2 ) equipped 
with the elastic metric that connects the cat-shaped figure 
to the dog-shaped figure. Original image published in [5]. 



A scale invariant version of the .fP-metric 



G c {h,k) 



l - [ (D s h,D s k)ds 

■c JS 1 



1301. There the authors derive 



4 j s 

has been studied in 
an explicit solution formula for the geodesic equation 
and calculate the sectional curvature. More general and 
higher order Sobolev metrics on plane curves have been 



(10) studied in [80j|90], and they have been applied to the 



field of active contours in 26.118 . Other Sobolev type 



metrics on curves that have been studied include a met- 
ric for which translations, scale changes and deforma- 
tions of the curve are orthogonal |116| and a ff 2 -type 
metric which has translations, scalings and rotations in 



its kernel 110 



For curves we can use arclength to identify each cl- 
ement C £ Bj,j(S ,M, 2 ) of shape space with a (up to 
rotation) unique parametrized curve c € Imm(5' 1 ,K 2 ). 
This observation has been used by Preston to induce 
a Riemannian metric on the shape space of unpara- 
metrized curves, via metrics on the space of arclength 
parametrized curves; see [103 104 . A similar approach 
has been chosen in (70l . 



6.2 Geodesic distance 

Sobolev-type metrics induce a point-separating geode- 
sic distance function on Bij(M, K rf ), if the order of the 
operator L is high enough. For the f/^-metric 

Gf(h,k)=[ ((Id+^)M)vol(<fs), (11) 

J M 

one can bound the length of a path by the area swept- 
out, similarly to the case of almost local metrics. 

Theorem 6.2 If the metric G L induced by the operator 



L is stronger than the H 1 -metric (11), i.e., 



G^(h,k)>CGf(h,k) 
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for some constant C > 0, then G L induces a point- 
separating geodesic distance function on the shape space 
B itf {M,R d ). 

A proof can be found in [11 | Thm 7.6]. An ingredient 
in the proof is the Lipschitz continuity of ^/Vol q . 

Theorem 6.3 The H 1 -metric satisfies 



1 



V^ol^-yVolQ 
In particular the map 
VWl; {B iif {M, K d ),distf^) 
is Lipschitz continuous. 



< ^distf* (Q U Q 2 ). 



90 



11 



Sec. 
Lem. 



A proof for plane curves can be found in 
4.7] and for higher dimensional surfaces in 
7.5]. 

The behavior of the geodesic distance o n th e space 
Imm/(M, M. d ) is unknown. Similar to Sec. 5.1 we can 
restrict the G L -metric to an orbit qoT)iff(M) and the 
induced metric on Diff(Af) will be a right-invariant 
Sobolev metric. Since Sobolev-type metrics of a suffi- 
ciently high order on the diffcomorphism group have 
point-separating geodesic distance functions, there is 
no a-priori obstacle for the distance dist^^ not to be 
point-separating. 

Open question. Under what conditions on the oper- 
ator L does the metric G L induce a point-separating 
geodesic distance function on Immj(M, M. d )7 



Theorem 6.4 (Thm. 6.5, (TT)) Let L be a smooth 
pseudo-differential operator, that is invariant under re- 
parametrizations, such that the adjoint Adj(VL) exists 



in the sense of (12). Then the geodesic equation for the 
G -metric on Imm(M, R d ) is given by: 

p = L q q t ® vol(q*g) 

pt = ^Adi(V L)(q u qt^ ~2Tq.(L q q t ,Vq t ) t (1 3 ) 
- H q (L q qt: qt)n q ) ® vo\(q*g) 



Note, that only the normal part of the adjoint 

Adj(VL)(gt,g t ) x = (Adj{V L)(q u q t ),n q )n q 

appears in the geodesic equation. The tangential part is 
determined by the reparametrization invariance of the 
operator L; see 11 Lem. 6.2]. 



Example 6.5 Consider the simple operator L = D s on 
the space Imm(S' 1 ,R 2 ) of plane curves. To emphasize 
the nonlinear dependence of L on the footpoint c we 
write it as L r h — rKrh' . The covariant derivative VL is 

|c'| 

simply the derivative of L with respect to the footpoint, 



V k Lh= —(h?,d)ti = (D s h,v c }D s h. 

for the operator L = D s . To compute its adjoint, we 
use the identity, obtained by partial integration, 



6.3 The geodesic equation 

The most concise way to write the geodesic equation 
on Imm(M, M. d ) for a general operator L involves its 
covariant derivative VL and adjoint Adj(VL). For a 
general operator L on Imm(M, R d ) we define the ad- 
joint Adj(VL) to be the adjoint of VkLh in the k vari- 
able with respect to the L 2 -metric, i.e., for all h,k,m € 
T q Imm(M, R d ) we have 

/ {V k Lh, m) vol(q*g) = 

= f (k,Adj(VL)(h,m))vol(q*g). (12) 

The existence and smoothness of the adjoint has to be 
checked for each metric by hand. This usually involves 
partial integration and even for simple operators like 
L = Id+(Z\ 9 )' the expressions for the adjoint quickly 
become unwieldy. 

Assuming the adjoint in the above sense exists, we 
can write the geodesic equation in the following form 
in terms of the momentum. 



/ (D s k,v c )(D s h,m) ds = 
Js 1 

= - (k,K c n c )(D s h,m) + (k,v c )D s (D s h,m) ds , 
Js 1 

which leads to 

Adj(VL)(/j,m) = -{D s h,m)n c n c - D s ({D s h,m))v c . 

The normal part Adj(VL)^, which is necessary for the 
geodesic equation is 

Adj(Vi) _L (/i, m) = —(D s h, m)n c n c . 

Note that while the full adjoint is a second order dif- 
ferential operator, the normal part has only order one. 
This reduction in order will be important for the well- 
posedness of the geodesic equation. 

To prove that geodesies on Bij(M,R d ) can be rep- 
resented by horizontal geodesies on Immj(M, IR'') we 
need the following lifting property. 
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Lemma 6.6 (Lem 6.8 and 6.9, (llj) Let L be a 

smooth pseudo-differential operator, that is invariant 
under reparametrizations, such that for each q, the op- 
erator L q is elliptic. Then the decomposition 

T Imm/ (M, R d ) = Hor L © Ver 

of tangent vectors into horizontal and vertical parts is 
a smooth operation. 

For any smooth path q(t) in Imm/ (M, R d ) there ex- 
ists a smooth path ip(t) in Diff (M) depending smoothly 
on q(t) such that the path q(t) = q(t) o ip(t) is horizon- 
tal, i.e., 

Gi {t) (d t q(t) ,Tq.X)=0, Vie X(M) . 

Thus any path in shape space can be lifted to a hor- 
izontal path of immersions. 



6.4 Well-posedness of the geodesic equation 

The well-posedness of the geodesic equation can be pro- 
ven under rather general assumptions on the operator. 

Assumptions. For fixed q G Imm(Af, R d ) the op- 
erator L q is an elliptic, pseudo-differential operator of 
order 21 and it is positive and symmetric with respect 
to the L 2 -metric. 

The operator L, the covariant derivative VL and 
the normal part of the adjoint Adj(VL) 1 - are smooth 
sections of the corresponding bundles. For fixed q they 

L q h, VL q (h,k), Adj(V L^k) 1 - 

are pseudo-differential operators of order 21 in h, k sep- 
arately. As mappings in the footpoint q they can be 
a composition of non-linear differential operators and 
linear pseudo-differential operators as long as the total 
order is less than 21. 

The operator L is reparametrization invariant in the 
sense of 

With these assumptions we have the following the- 
orem from Thm. 6.6]. A similar theorem has been 
proven for plane curves in |90| Thm 4.3]. 

Theorem 6.7 Let the operator L satisfy the above as- 
sumptions with I > 1 and let k > dlm ( M *> + 21 + 1. Then 
the geodesic spray of the G L -metric is smooth on the 
Sobolev manifold lm.m k (M,R d ) of H k -immersions. 
In particular the initial value problem for the geode- 



sic equation ( 13 I has unique solutions 



Remark 6.8 For smooth initial conditions q(0, •), q t (0, •) 
in Tlmm k (M,R d ) we can apply the above theorem for 
different k and obtain solutions in each Sobolev comple- 
tion Imm fe (M, R d ). It can be shown that the maximal 
interval of existence is independent of the Sobolev order 
k and thus the solution of the geodesic equation itself is 
in fact smooth. Therefore the above theorem continues 
to hold, if Imm fe (M,K d ) is replaced by Imm(A/, R d ). 

Remark 6.9 Due to the correspondence of horizontal 
geodesies on Imm/(M, R d ) and geodesies on on shape 
space Bij(M, R d ) the above well-posedness theorem 
implies in particular the well-posedness of the geode- 
sic problem on B hf (M, R d ). 

Example 6.10 The assumptions of this theorem might 
look very abstract at first. The simplest operator ful- 
filling them is 

L q =ld+A a 

or any power of the Laplacian, L q = Id+(Z\ 9 )'. We can 
also introduce non-constant coefficients, for example 

L q = f 1 (H q ,K q ) + f 2 (H q ,K q )A\ 

as long as the operator remains elliptic, symmetric and 
positive. To check symmetry and positivity it is some- 
times easier to start with the metric, for example 



G q {h,k)= / gi (H q ,K q )(h,k) + 

JM 

d 

+ g 2 {H q , K q ) ]T g(V 9 h, VW) vol(q*g) . 

i=l 

From this we can compute the operator L q via partial 
integration 

(L q hY - 9l (H q ,K q )h* - div 9 (g 2 (H q ,K q )W 9 h l ) . 

We see that this operator only satisfies the assumptions 
if <?2 is the constant function, because for fixed q the 
operator L q h has order 2 in h, but it would depend on 
the third derivatives on q. 



6.5 Conserved quantities 

If the operator L is equivariant with respect to repa- 
rametrizations, the G L -metric will be invariant under 
the action of Diff(M). By Noether's theorem the repa- 
rametrization momentum is constant along each geode- 



t^q(t,-) e C 00 ((-e,e),Imm fe (Af,]R d )), 

for small times and the solution depends smoothly on 
the initial conditions q(0, •), q t (0,-) in rimm"(M, R d ). 



sic, c.f. Sect. 5.3 This means that for each X 6 £(M) 
we have 



/ (L q q t ,Tq.X) vo\(q* g) = const. 
JM 
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If L is additionally invariant under the action of the 
Euclidean motion group M. d x SO(d) then so is the G L - 
metric and the linear and angular momenta are con- 
stant along geodesies. These are given by 



(L q q t ) vol(q*g) G 



M 



1 qA(L q q t )val(q*g) G /\" 

M 



so(dy 



If the operator L satisfies the scaling property 



X 



- dim(M)-2 7 



q G Imm(M, 



,AG 



then the induced metric G L is scale invariant. In this 
case the scaling momentum is conserved along geodesies 
as well. It is given by: 



(q, L q q t ) vol (3) G 



M 



See Sect. 5.3 for a more detailed explanation of the 



meaning of this quantities. 



6.6 Completeness 

Concerning geodesic completeness it is possible to de- 
rive a result similar to Thm. 15.71 The set of concentric 
spheres with a common center is again a totally geode- 
sic submanifold and we can look for conditions, when it 
is possible to shrink spheres to a point with a geodesic 
of finite length. 



Theorem 6.11 (Lem. 9.5, (TTl) If L = If+(A 



and I < 



i(Af) 



2 



+ 1. then the spaces Imm(5 d 1 ,R ) 
d ) are not geodesically complete with 



and Bij(S' 
respect to the G^ -metric. 

For other choices of M scalings will in general not 
be geodesic, but under the same condition an immer- 
sion can be scaled down to a point with finite energy. 
Under what conditions are spaces become geodesically 
complete is unknown. We do however suspect that sim- 



ilarly as Thm. 7.5 for the diffeomorphism group, a dif- 
ferential operator of high enough order will induce a 
geodesically complete metric. 

The metric completeness of Bij(S 1 , K 2 ) is known 
for the Sobolev metrics 



1 



G^(h,k) = I ^{h,k)+£ 2l {Dlh,D{k)As, 

with j = 1,2. For the metric of order 1 we have the 
following theorem. 

Theorem 6.12 (Thms. 26 and 27, (8(jj) The metric 
completion of the space (Bij, dist H ) is B l j lp (S 1 , R 2 ), 
the space of all rectifiable curves. 



See Sect. 5.4 for details about B^ P (S 1 ,M 2 ). There 



is a similar result for the metric of second order. 



Theorem 6.13 (Thm. 29, [80]) The completion of 

H 2 \ 

the metric space [Bij, dist ) is the set of all those 
rectifiable curves that admit curvature k c as a measur- 
able function and J sl n 2 ds < oo. 



6.7 Curvature 

Apart from some results on first and second order met- 
rics on the space of plane curves, very little is known 
about the curvature of Sobolev-type metrics on either 



Imm/(M,R d ) or B itf (M,R d ) 



The family ( 10 ) of if 1 -type metrics on the space 
lmm([0, 27r], E 2 )/trans of open curves modulo transla- 
tions is isometric to an open subset of a vector space 
and therefore flat; see [5], It then follows from O'Neil's 
formula that the quotient space Bit(M, M. d ) has non- 
negative sectional curvature. 

For the scale-invariant H^-type metric 



G c {h,k) = — / (D s h,D s k)ds, 
H-c Js 1 

the sectional curvature has been computed explicitly on 
the space B i j(S 1 ,M 2 )/ (aim), which is the quotient of 
Bjj(S' 1 ,]R 2 ) by similarity transformations — transla- 
tions, rotations and scalings — in 130 . The sectional 



curvature is again non-negative and upper bounds of 
the following form can be derived. 



Theorem 6.14 (Sect. 5.8, [130]) Take a curve c G 
Imm(S' 1 ,IR 2 ) and let hi,h 2 G T c Imm(5 1 , R 2 ) be two 
orthonormal tangent vectors. Then the sectional curva- 
ture atC — 7r(c) of the plane spanned by the projections 
of '/ii,/i2 in the space B^t^S , ]R 2 )/(sim) is bounded by 

< k c (P(h 1 ,h 2 )) <2 + A 1 (k c )+ 

+ A 2 (K c )\\D s h 2 ■ n||oo + A 3 {K c )\\D s (D s h 2 ■ n)^ . 

where Ai : C°°(5' 1 ,K) — > K are functions of K that 
are invariant under reparametrizations and similarity 
transformations. 



Explicit formulas of Ai(n) can be found in 130 
This is a bound on the sectional curvature, that depends 
on the first two derivatives of h 2 and is independent of 
hi . Thus geodesies have at least a small interval before 
they meet other geodesies. 



A similar bound has been derived in 110 for a sec- 



ond order metric on the space of plane curves. 
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7 Diffeomorphism groups 

In the context of shape spaces diffeomorphism groups 
arise two-fold: 

— The shape space B i f(M,W l ) of immersed subman- 
ifolds is the quotient 

B itf (M,R d ) = Imm / (M,R d )/Diff(M) 

of the space of immersions by the reparametrization 
group Diff(M). 

— By fixing an embedding go € Emb(M, M. d ) we have 
the map 

Diff c (M d ) Emb(M, R d ), ^p 9o 

A right-invariant Riemannian metric on Diff c (IR d ) 
induces a Riemannian metric on Emb(M, R d ) such 
that the above map is a Riemannian submersion. 
See Sect. [H for this construction. 

These are the two main applications of the diffeomor- 
phism group discussed in this paper. Thus we mainly 
will treat the group 

Diff(M) = {ipeC°°(M,M) : tp bij., <^ _1 € C°° (M, M)} 

of smooth diffeomorphisms of a closed manifold M and 
Diff (R. d ) , groups of diffeomorphisms of R d . On WL d we 
will consider diffeomorphism groups with the following 
decay conditions towards infinity 



Diff c (R d ) 
Diff 5(1^) 
Diff H oo(]R d ) 



{if : supp(<p — Id) compact} 
: (ip-M)eS(R d )} 



Here H°°(M. d ) denotes the intersection of all Sobolev 
space H k (R d ), k e N> and S(R d ) denotes the space 
of all rapidly decreasing functions. They are smooth 
regular Lie groups. Their Lie algebras are the spaces 
X c (K d ),X 5 (K d ) and X H °°{R d ) of compactly supported, 
rapidly decreasing and Sobolev vector fields, respec- 



tively. See 86 and 91 for details 



On general non-compact manifolds N one can also 
consider the group of compactly supported diffeomor- 
phisms like on R d ; see 73 Sect. 43]. 



7.1 Right-invariant Riemannian metrics 

A right-invariant metric on Diff c (K <i ) is determined via 

G V {X V ,Y V ) = (I.o^^y.o^i, 

X V ,Y V e T v DiS c (R d ) , (14) 



by an inner product (•, -)l on the space X c (18i d ) of vec- 
tor fields. We assume that the inner product is defined 
via a symmetric, positive definite, pseudo-differential 
operator L : X(R d ) -> X(R d ) by 



(X,Y)l = 



LX -Ydx. 



Examples of such inner products include 

— The L 2 -metric with L = Id, 

(X,Y)& = / (X,Y)dx. 

— The Sobolev-type metrics of order s with s > 0, 
(X,Y) H .= [ (l + |e| 2 ) S (A(0?(0)de, 



with X(0 = (2Tr)- d / 2 f Rd e-' ll - x ^X{x) dx being the 
Fourier transform. Note that for s£N these metrics 
can be written as 

(X,Y) H s = J^((ld-A) S X,Y) dx, 



i.e., L = (Id- A) 
The fami 

(X,Y) a>h 



The family of a-&-c-metrics, introduced in 66 
/ a(X,Y) 

JR d 



6divXdivF+ 
+ c(VX,VY)dx. (15) 



Recall that A = —A 9 denotes the usual Laplacian on 
M. d , which is the negative of the geometric Laplacian; 
see Sec l2.ll In dimension d — 1 the second and the third 
term coincide and the metric simplifies to the family of 
a-b metrics 

(X,Y) a:b = [ aX ■ Y + bX'Y 1 dx . 
Jr 

On manifolds other than M. d , one can use the in- 
trinsic differential operators to define inner products 
on X(M), which are then extended to right-invariant 



Riemannian metrics on Diff(M) via (14). For exam- 



ple, when (M, g) is a Riemannian manifold the Sobolev- 
metrics of integral order can be defined using the Lapla- 
cian A 9 via 

(X,Y) Hk = [ g((Id+A 9 ) k X,Y)vol a . 

J M 

Similarly the family of a-b-c metrics have an intrinsic 
representation given by 



(X,Y) aAc = / ag{X,Y) + bg{5X b ,5Y*)+ 

JM 

+ cg(dX b ,dY b )vo\ 9 
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More general Sobolev spaces H S (M) with s / N and 
the corresponding norms can be introduced using parti- 
tions of unity and Riemannian exponential coordinates. 
See the books 119] and [43] for the theory of function 
spaces, including Sobolev spaces of fractional order, on 
manifolds. 

Remark 7. 1 An alternative approach to induce a metric 
on the diffeomorphism group is to use a reproducing 
kernel Hilbert space H of vector fields, with X(R d ) C U 
and consider the restriction of the inner product on H 
to This approach is described in Sect. 8.2 



7.2 Geodesic Equation 

The geodesic equation on any Lie group Q with a right- 
invariant metric is given as follows. A curve g(t) 6 Q 
is a geodesic if the right logarithmic derivative u(t) = 
<9 t g(i)g(i) _1 satisfies 

dfU = — ad^ u , 

where ad T is the transpose of ad with respect to the 
given inner product •) on the Lie algebra, i.e., 

7(ad^, v, w) = y(v, ad M w) 

On DifF(R d ) with a metric given via an operator L 
we can write the equation as a PDE in terms of the 
momentum m = Lu, 

dtm + (u- V)m + m div u + Du T .m = 0, m = Lu , (16) 

and dt^p = uoip. For different choices of L one can 
obtain the following PDEs as geodesic equations. 

The L 2 -metric with Lu — u in one dimension has as 
geodesic equation Burgers' equation, 

u t + 3uu x = . 

This equation was used as a model equation for turbu- 
lence in 23 . 

The H 1 -metric with Lu = u — u xx in one dimension 
has as geodesic equation the Camassa-Holm equation 
221, 



u t - u xxt + 3uu x - 2u x u xx 



= 0. 



It describes the propagation of shallow water waves on 
the free surface under gravity. It is a completely inte- 
grable equation and possesses a bihamiltonian struc- 
ture, that gives rise to an infinite number of conserva- 
tion laws. 

The homogeneous if ^metric on DiS(S 1 ) with the 
operator Lu — —u xx has as geodesic equation the peri- 
odic Hunter-Saxton equation 



Since the IP-metric has a kernel, we have to consider it 
as a metric on the right coset space Rot(S ,1 )\ Diff (S* 1 ) 
of diffeomorphisms modulo rigid rotations. The Hunter- 
Saxton equation was proposed as a model for the prop- 
agation of orientation waves in nematic liquid crystals 
in [57]. Its geodesic nature was discovered in 68 . It 
is also a completely integrable, bihamiltonian equation 
with an infinite number of conservation laws 58 . As a 



Riemannian manifold (Diff (S 1 ), H 1 ) is isometric to an 
open subset of a sphere and as such has positive con- 
stant curvature 76 . It was shown recently in j9j, that 
a related result also holds for the non-periodic Hunter- 
Saxton equation, which is the geodesic equation for the 
7f 1 -metric on a certain extension of Diff c (M). 

Between the Hunter-Saxton and the Camassa-Holm 
equation lies the fi-Hunter-Saxton equation, 

u xxt - 2fi(u)u x + 2u x u xx + uu xxx = , 

which is the geodesic equation on the circle with re- 
spect to the ^fi^-metric defined by the operator Lu = 
fj.(u) — u xx , with n(u) = 5^ J s i udx being the mean. It 
was introduced in [65^ as a non degenerate metric on 
Diff(S' 1 ), such that the projection 

(Diff(S' 1 ),^ij 1 ) -> (Rot(S* 1 )\Diff(S ,1 ),ii 1 ) 

is a Riemannian submersion. It is also a completely in- 
tegrable, bihamiltonian equation. 

The geodesic equation for the homogeneous H 1 / 2 - 
metric is the modified Constantin-Lax-Majda (mCLM) 
equation, 



m t + um x + 2u x m = 0, 



Hu x 



The mCLM equation is part of a family of one dimen- 
sional models for the vorticity equation [36 39 101 
Its geodesic nature was recognized in 



123 . As for the 



Hunter-Saxton equation we have to regard the mCLM 
equation on the coset space Rot(S' 1 )\ Diff (S 1 ). 

In the context of hydrodynamics a closely related 
space is the Virasoro-Bott group 

Vir(S* 1 ) =Mx c Diff(5 1 ), 



with the group operations 
a)\p) \a + /3 + c(ip, $)) ' 



V" 1 
—a 



^xxt H~ ^^x^xx ~t~ tilt's. 



0. 



for tp, tp e Diff(5 1 ), and a,P € R. The Virasoro-Bott 
group a central extension of Diff(S' 1 ) with respect to 
the Bott-cocycle: 

c : Diff(S' 1 ) x Diff(5 1 ) ->■ R 

c(</5,V) = 2 flogfy o^)dlogV>', 
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Fig. 7 Some geodesic equations on diffeomorphism groups, 
that are relevant in mathematical physics. 



and it is the unqiue non-trivial central extension of 
Diff (S 1 ). For a detailed exposition of the Virasoro-Bott 
group see the book of Guieu and Roger [54]. It was 
found in 102 108 that the geodesic equation of the 



right invariant L 2 -metric on the Virasoro-Bott group is 
the Korteweg-de Vries equation 



Ut + 3u x ii + au x 



0, a G 



Similarly the Camassa Holm equation with dispersion 

Uf — Uixxt ~ t - QuU% — ^UxUxx — ^^xxx ~ t - 2/KrUx 

was recognized to be the geodesic equation on the Vira- 



soro-Bott group with respect to the i? 1 -metric in 96 . 



The geodesic equation ( 16 1 can be rewritten as an 



integral in Lagrangian coordinates. For a metric given 
by a differential operator, let K(x, y) be its Green's 
function. We assume that the initial momentum mo 
is a vector-valued distribution, whose components are 
finite measures. The initial velocity can be obtained 
from mo via Uq(x) = J Rd K(x, .)m (.) and conversely 
Luq <8> dx. The Lagrangian form of ( 16 ) is 



m 



d t (f(t,x) 



K(cp(t, x),ip(t, .)) ip(t)*m (.) . 



7.3 Well-posedness of the geodesic equation 



One possible method to prove the well-posedness of the 
geodesic equations is to extend the group and the metric 
to the Sobolev-completion 

Diff 9 (M) = {if G H q (M, M) : ip\A]., v?" 1 G H q (M, M)} 

which is a Hilbert manifold and a topological group 
for q > dim(Af)/2 + 1. It is however not a Lie group 
any more, since the right-multiplication is only con- 
tinuous but not smooth. Nevertheless it is possible to 
show that the geodesic spray of various metrics on the 
Sobolev-completion is smooth for q large enough and 
then an application of the theorem of Picard-Lindelof 



for ODEs shows the existence and smoothness of the 
exponential map. This method was first applied in 42 
for the L 2 -metric on the group of volume-preserving 
diffcomorphisms to show the existence of solutions for 
Euler's equations, which model inviscid, incompressible 
fluid flows. On the full diffeomorphism group the follow- 
ing well-posedness results can be obtained via the same 
method. 

Theorem 7.2 (Thm. 3.3, (48)) Let (M,g) be a com- 
pact Riemannian manifold without boundary. The geo- 
desic spray of the H 1 -metric 

(u,v) = / g(u,v) + g(Vu, Vv) vol 9 

J M 



is smooth as a map TDiff 9 (M) T 2 Diff 9 (M) for q > 

dim(M) 
2 



1. 



The (higher- dimensional) Camassa-Holm equation 
with initial condition uq G X q (M,Al) admits a unique 
solution u(t) for small times and we have 

t^u(t) e C {(~e,e),X q (M))DC 1 ((~e,e),X q - 1 (M)), 

and the map uq G X q (M) H> u(t) G X q (M) is continu- 
ous. 

This result holds more generally also for manifolds 
with boundary with either Dirichlet, Navier or mixed 
boundary conditions. See [48] for more details. For the 
one-dimensional case the smoothness of the geodesic 



spray was noted already in 71 

For the circle M — S 1 we have the stronger result 
that the geodesic sprays for Sobolev metrics H s are 
smooth for s > ^. 

Theorem 7.3 (Cor. 4.2, |44|) The geodesic spray of 
the H s -metric 



( 17 ) (it, v) = ^(1 + n 2 ) s u{n)v{n) 



on the diffeomorphism group Diff 9 (5 1 ) of the circle is 
smooth for s > | and q > 2s + | . Thus the geodesic 
equation is, similarly to Thm. \ 7.2\ locally well-posed. 

The case of Sobolev metrics of integer order, which 
include s th e periodic Camassa-Holm equation, was pro- 
ven in |35| . For the homogeneous ij^-metric this re- 
sult was proven in [45] and the estimates were then ex- 
tended to cover general metrics given via Fourier mul- 



tipliers in 44 



As a consequence of the well-posedness result for So- 
bolev metrics on Imm(M, M d ) it has been shown that 
the Lagrangian form of the geodesic equation is lo- 
cally well-posedness for higher order Sobolev metrics 
on Diff(M). 
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Theorem 7.4 (Thm. 10, [11]) Let (M,g) be a com- 
pact Riemannian manifold without boundary and let 
k > 1 with keN. Forq> dim 2 (M) + 2k + 1 the geodesic 
spray of the H k -metric is smooth as a map on Diff 9 (M) 
and the geodesic equation has unique local solutions on 
Diff 9 (M). 

If the metric is strong enough, it is possible to show 
the long-time existence of solutions. 

Theorem 7.5 // the Green's function K of the opera- 
tor L inducing the metric is a C 1 -function, then for any 
vector-valued distribution mo, whose components are fi- 
nite signed measures, equation (171 with ip(0,x) = x 



can be solved for all time and the solution is a map 
*i->¥>(*.0 G C' 1 (M,C' 1 (M d ,M d )). 



This result is implicit in the work 121 , an explicit 



proof can be found in 100 . See also 129 . 



Remark 7.6 This method of proving well-posedness is 
not universally applicable as not all geodesic sprays are 
smooth. For example the spray induced by the right- 
invariant L 2 -mctric on Diff(S' 1 ) is not smooth. More 



precisely in 34 it is shown that the exponential map 



is not a C 1 -map from a neighborhood of Tjd Diff 9 (5' 1 ) 
to Diff 9 (5 1 ) for any q > 2. Nevertheless the geodesic 
equation, which is Burgers' equation in this case, has 
solutions for small time 

t^u(t)eC ((-e,e),H q {S 1 ))nC 1 ((~e,e),H q - 1 {S 1 )), 



when uq £ H q (S 1 ); see 63 . A similar statement holds 



for the KdV-equation, which is the geodesic equation 
on the Virasoro-Bott group with respect to the right- 



invariant L 2 -metric; see 33 . 



7.4 Geodesic distance 

It was shown in [88] that the geodesic distance on the 
group Diff C (M) vanishes for the L 2 -metric and is posi- 
tive for the iJ 1 -metric. This naturally raises the ques- 
tion, what happens for the ff s -metric with < s < 1. 
For M — S 1 a complete answer is provided in [6], 
whereas for more general manifolds N a partial answer 
was given in the articles [6j[8j. 

Theorem 7.7 (Thm. 3.1, (6[[8)) The geodesic dis- 
tance on Diff c (M) induced by the Sobolev-type metric 
of order s vanishes 

- for s < \, 

— for s = 5, when M — S 1 x C with C compact. 
The geodesic distance is positive 



— for s > 1, 

— for s > |, when dim(Af) = 1. 

Remark 7.8 By taking C = {*} to be the zero dimen- 
sional manifold, the above theorem provides a complete 
answer for M = S : the geodesic distance vanishes if 
and only if s < | . 

Remark 7.9 The _ff ^-metric on Diff(S' 1 ) is the only 
known example, where the geodesic spray is smooth on 
the Sobolev-completions Diff 9 (S' 1 ) for all q > | and the 
geodesic distance vanishes at the same time. 

It is shown in [44] that for q > | the exponential 
map is a local diffeomorphism exp : U C H q (M.) — > 
Diff 9 (S 1 ) . In particular we have the inequality 



Len' 



'(¥>)> Ijexp- 1 ^!)))!^ 



for all paths ip : [0, 1] exp(B|(0)) with tp(0) = Id. In 
other words we have a lower bound on the length for all 
paths, that remain H g -c\ose to Id. This does not how- 
ever imply anything about the geodesic distance, since 
a path can have small iJ 1 / 2 -length or equivalently re- 
main H 1 / 2 -close to Id, but leave the ^-neighborhood. 

Open question. For a Sobolev metric of order s the 
behavior of the geodesic distance on Diff c (M) remains 
open for 

- 1 < s < 1 and dim(iV) > 2. 

- s = \ and N ^ S 1 x M, with M compact. 

Extrapolating from the known cases, we conjecture the 
following result: The geodesic distance induced by the 
Sobolev-type metric of order s onDiS c (N) vanishes for 
s < o and is non- degenerated for s > 5- 

A main ingredient for the vanishing result is the fol- 
lowing property of the geodesic distance on Diff c (iV): 

Lemma 7.10 If the geodesic distance on Diff c (iV) van- 
ishes for one pair ip, if) £ Diff c (iV) with ip ^ tj), then the 
geodesic distance already vanishes identically on all of 
Diff c (A0. 

This lemma follows from the fact that the set 



A = 



{ip : dist ffS (Id, ip) = 0} 



is a normal subgroup of Diff c (TV) for all s and because 
Diff C (N) is a simple group. Thus, if A contains any 
element apart from Id it has to be the whole group. 

Remark 7.11 We can also consider the geodesic dis- 
tance on the Virasoro-Bott group, which is the one- 
dimensional central extension of Diff(S' 1 ). There the 
geodesic distance vanishes for s = 0, i.e., for the L 2 - 
metric. For s > | the geodesic distance cannot vanish 
identically. Whether it is point-separating definite is not 
known. 
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Open question. For a Sobolev metric of order s the 
behavior of the geodesic distance on the Virasoro-Bott 
group remains open for < s < 1. 

One way to define geodesies is to fix two diffcomor- 
phisms ifo, ipi and to consider the set 

B = {<p(t) : <p(0) = tpo, tp(l) = <M 

of all paths joining them. Geodesies then correspond to 
critical points of the energy or equivalently the length 
functional restricted to the set B. Vanishing of the geo- 
desic distance implies that these functionals have no 
global minima. The following theorem shows that for 
the L 2 -metric there are no local minima either. 

Theorem 7.12 (Thm 3.1, (2l)) Let ip(t,x) with t e 
[0, T] be a path in Diff c (M). Let U be a neighborhood of 
ip in the space C£°([0, T] xR). Then there exists a path 
ip € U with the same endpoints as ip and 

E($) < E(ip), 

where E(.) is the energy w.r.t. the right-invariant L 2 - 
metric. 

In the article [21] the result is proven for Diffs(R), but 
it holds for DifLflR) as well. 



7.5 Completeness 

As a corollary of Thm. [73] we obtain the result that the 
diffcomorphism equipped with a metric of high enough 
order is geodesically complete: 

Theorem 7.13 Let (M,g) be a compact Riemannian 
manifold and let G s be the Sobolev metric of order s. 



For s > 



dim(Af)+3 



the space (Diff(M),G s ) is geodesi- 



cally complete. 



This result is based on the observation, that for s > 
dim(M)+3 kernel of the metric G s is a C 1 -function. 



7.6 Curvature 

Denote by •) the inner product on the Lie algebra g 
of any Lie group Q and let u, v € q be orthonormal vec- 
tors. Then the sectional curvature of the plane P(u, v) 
in Q with respect to the right-invariant metric induced 
by 7 is given by 



k(P(u, v)) 



-\\ a,dv u + &dl v\\ 2 
1 



7(ad^ u, adl v) 



where ad is the transpose of ad with respect to the 
given inner product 7 inducing the right invariant met- 
ric. 

For general Sobolev metrics there are no results on 
curvature available, but for the family of a-6-c-metrics 
( 15 ) on the d-dimensional torus T d , it was shown in 66 



that the curvature assumes both signs. 



|| ad„ v\\ + -7(ad u v, ad v u - ad u v), 



Theorem 7.14 (Thm. 7.1, [66]) If d > 2 and at 

least two of the parameters a, b, c are non-zero, then 
the sectional curvature of the a-b-c-metric on Diff(T d ) 
assumes both signs. 

In dimension one we have the same behavior for the 
family of a-b metrics. 

Theorem 7.15 (Sect. 6, [66]) If d = 1 and both pa- 
rameters a, b are non-zero, then the sectional curvature 
of the a-b-metric on Diff (S 1 ) assumesboth signs. 

There are two special cases, where the sign of the cur- 
vature is constant. The first is the L 2 -metric (b = 0) in 
one dimension. 

Theorem 7.16 (Sect. 5.4, (88)) If d = 1 and b = 

then the sectional curvature of the plane spanned by 
two orthonormal vector fields u, v £ X(S 1 ) for the a- 
b metric on Diff(S' 1 ) is given by 

KP M )=( (uv'-vu'fte. 
Js 1 

In particular the sectional curvature is non-negative. 
This does not generalize to higher dimensions. 

Theorem 7.17 (Prop. 7.2, [66]) If d > 2 and b = 

c = then then the sectional curvature of the a-b-c 
metric on Diff(T d ) assumes both signs. 

The sectional curvature of the L 2 -metric has been cal- 
culated for an arbitrary Riemannian manifold N. The 
expression for sectional curvature is the sum of a non- 
negative term and a term whose sign is indefinite. Al- 
though we conjecture that the statement of Thm. |7.17| 
extends to arbitrary manifolds N, this has not been 
proven yet. 

The second special case is the homogeneous H 1 - 
metric with a = c = 0ford>2 and a = for d = 1. 
The metric is degenerate on Diff(M), but it induces the 
Fisher-Rao metric on the space Diff (M)/ Diff ^(M) of 
densities. Remarkably the induced metric has constant 
sectional curvature. 

Theorem 7.18 (Cor. 3.2, (66)) Let (M,g) be a com- 
pact Riemannian manifold. Then the homogeneous H 1 - 
metric 

(u, v)jji = / div(u) div(w) vol(g) 
Jm 
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on Diff(Af)/ Diff^(M) has constant positive sectional 
curvature 

This result is based on the observation, that the H 1 
metric on DirF(M)/Diff Al (M) is isometric to a sphere 
in the Hilbert space L 2 (M, vol(g)). For M = S 1 this 
result has been proven already in [76]. Recently it has 
been shown that the H 1 -metric on a certain extension 
of Diff c (R) is a flat space in the sense of Riemannian 
geometry; see [9]. 

8 Metrics on shape space induced by Diff(M Tl ) 

In this section we will consider Riemannian metrics on 
B e (M,R d ), the space of embedded type M subman- 
ifolds that are induced by the left action of Diff(]R' i ). 
The group Diff(R d ) can be any of the groups Diff c (R d ), 
Diff 5 (K d ) or Diff ff oo(M d ) described in Sect.0 We also 
relax the assumption on the dimension of M and only 
require dim(M) < d. The action is given by 

Diff(M d ) x B e (M,R d ) 3 (ip,Q) h-> tp(Q) e B e {M,R d ) . 

This action is in general not transitive, but its orbits 
are open subsets of B e (M, Mr). If we restrict the action 
to Diffo , the connected component of the identity, 
the orbits are the connected components of B e (M, R d ). 
For Q e B e (M,R d ) the isotropy group 

Diff (R d ) Q = {<p : tp(Q) = Q} , 

consists of all diffeomorphisms that map Q to itself. 
Thus each orbit Orb(Q) = Diff can be identified 

with the quotient 

B e (M,R d ) D Orb(Q) = Diff(M <i )/Diff (R d ) Q . 

Let us take a step backwards and remember that an- 
other way to represent B e (M, R d ) was as the quotient 

B e (M,R d ) ^ Emb(M,K d )/Diff(M) . 

The diffeomorphism group Diff(]R d ) also acts on the 
space Emb(M, R d ) of embeddings, that is parametrized 
submanifolds with the action 

Diff (R d ) x Emb(M, R d ) 3 (cp, q) i-» <p o q e Emb(M, R d ). 

This action is generally not transitive either, but has 
open orbits as before. For fixed q G Emb(M, R d ), the 
isotropy group 

Diff(R d ), = {if : p\ q(M) =ld} , 



consists of all diffeomorphisms that fix the image q(M) 
pointwise. Note the subtle difference between the two 
groups Diff(R d ) g and Diff(R%, when Q = q(M). The 
former consists of diffeomorphisms that fix q{M) point- 
wise, while elements of the latter only fix q(M) as a 
set. As before we can identify each orbit Orb(q) = 
Diff(K rf ).g with the set 

Emb(A/,M d ) D Orb(g) Diff(R d )/Diff(R d )q . 

The isotropy groups are subgroups of each other 

Diff (R d ) q < Diff(M% < Diff (R d ) , 

with Diff(lR d ) g being a normal subgroup of Diff(R rf )Q. 
Their quotient can be identified with 

Diff {R d ) Q / Diff (R d ) q ~ Diff (M) . 

Now we have the two-step process, 

Diff(R d ) ->■ 

->■ Diff {R d )/ Diff {R d ) q S Orb(g) C Emb(M,M d ) 

-> Emb(M,M d )/Diff(Af) = B e (M,R d ) . 

In particular the open subset Orb(Q) of B e (M,R d ) can 
be represented as any of the quotients 

Orb(Q) £ Orb(g)/Diff(M) S 

~ Diff(K d )/Diff(IR d )„ / 

/ Diff (R%/ Diff (R d ) q ~ 

= Diff(E d )/Diff(M' i )Q . 

Let a right-invariant Riemannian metric G D be 
given on Diff(R d ). Then we can define a metric on 
Emb(M, R d ) in the following way. Fix q € Emb(M, R d ) 
and let q = tp o qo be an element in the orbit of qo . Then 
the norm of a tangent vector h € T q Emb(M, R d ) is de- 
fined by 

G^ h (h,h)= inf G™(X V ,X V ), 

X v oq = h 

with X v e T v Diff (R d ). If we define 7r go to be the pro- 
jection 

n qo : Diff(K d ) ->• Emb(M,M d ), n qo {p) = poq , 
then 

/;, = X v o q = T v TT qo .X v , 

and the equation defining the metric is the relation be- 
tween two metrics that are connected by a Riemannian 
submersion. Because G Dlff is right-invariant and the 
group action is associative we can rewrite the defining 
equation as 

G^ h (h,h) = inf G$ S (X,X), 
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with X e T M Diff(R d ). Thus we see that the metric 
G Emb is well-defined and does not depend on the choice 
of q$ . While ir qo is a Riemannian submersion this is an 
example, where the horizontal bundle exists only in a 
suitable Sobolev-completion; see Sect. |2.2| In Sect. |8.2| 
we will take care of this by defining the metric via a 
reproducing kernel Hilbcrt space %. 

This construction yields a Riemannian metric on the 
space Emb(-M, Mr) that is invariant under reparametri- 
zations, because the left-action by Diff(R d ) commutes 
with the right-action by Diff(M): 



G%%{ho<p,ho<p)= inf G^{X,X) = 

X o q o (p—fi o ip 

= inf G™ s (X,X) = G^ h (h,h). 

X o q=h ' 

The metric G Emb then projects to a Riemannian metric 
on B e (M,R d ) as explained in Sect. |2~2 



8.1 Pattern theory 

This section is closely related to ideas in Grenander's 
pattern theory 52|53|99 . The principle underlying pat- 
tern theory is to explain changes of shape by a defor- 
mation group acting on the shape. In our case shapes 
are elements of either Emb(M,K d ) or B e (M,R d ) and 
the deformation group is the group Diff(R d ). 

There is a lot of flexibility in the choice of the group 
and the space it acts upon. If M is a finite set of n 
points, then Emb(M,M d ) C (R d ) n is the set of land- 
marks. We have inclusion instead of equality because 
landmarks have to be distinct points. We will return to 
this space in Sec. [9| 

An important example is when the shape space is 
the space of volumetric grey-scale images modeled as 
functions in C°°(IR d ,IR) and the deformation group is 
Diff (R. d ) . The action is given by 



Diff(M d ) x C°°(I 



13(^,1) h> Io^ 1 G C°°(I 



This action is far from being transitive. Thus it is not 
possible to rigorously define a Riemannian metric on 



that is induces by Diff(K ). Nevertheless 



the idea of images being deformed by diffeomorphisms 
led to the image registration method known as LD- 
DMM 



8.2 Defining the metric on Diff (R d ) 

Following the presentation in [84] we assume that the 
inner product on X(R d ) is given in the following way: 
let (H, (•, be a Hilbert space of vector fields, such 
that the canonical inclusions in the following diagram 



X(I 



n ^ c£(i 



with k > are bounded linear mappings. We shall 
also assume that the Lie algebra X(M d ) of Diff(M d ) 
is dense in H. Here C^(Mr) is the space of all glob- 
ally bounded C fc -vector fields with the norm ||X||fc j0O = 

Under these assumptions % is a reproducing kernel 
Hilbert space, i.e., for all x, a G R d the directional point- 
evaluation ev° : % — > K defined as ev°(w) = (u(x),a) 
is a continuous linear functional on H. See [2 106] for 
a detailed treatment. The relation 

(u,K(.,x)a) n = (u(x),a) 

defines a matrix- valued function K : R d x R d — > R dxd , 
called the kernel of H. It satisfies the two properties 

- K G C%(R d ,R d ) and 

- K(y,x) = K(x,y) T . 

Associated to H we have the canonical isomorphism 
L : % —> H* . Note that the kernel satisfies K (y, x)a = 
L _1 (ev°)(y); this relation is even more general: The 
space 9Jt fc (lR d , of vector- valued distributions, whose 
components are fc-th derivatives of finite signed mea- 
sures is a subspace of the dual space C£(R d ,R d )* and 
the operator 

K : m k {R d ,R d ) -> Cj;(R d ), m^f K(.,x)m(x) 



coincides with L . This is represented in the diagram 



L 



-m k {R d ) 



c£(R d ,R d y 



K 



dy 
X(R d ) 



Here X(M d )* denotes the space of vector-valued distri- 
butions dual to X(M d ), depending on the decay condi- 
tions chosen. The inner product on X(K d ) is the restric- 
tion of the inner product on H, 



16||92]|93||120]. It is being applied in computa- (X, Y) H = (LX, Y) dx . 

;™„„ u„: jidt „„j r^T^ „„„ JR d 



tional anatomy with images being MRI and CT scans 
to study the connections between anatomical shape and 



physiological function. See 22 for an overview of image 



registration within the LDDMM framework. 



where the expression on the right hand side is a sugges- 
tive way to denote the pairing (LX, Y)x*xx between a 
distribution and a vector field. 
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Example 8.1 Let H = H k (R d ) be the Sobolev space of 
order k > | with the inner product 

(X,Y) H u= [ ((Id-A) k X,Y)dx. 

Then by the Sobolev embedding theorem we have 

H k (R d ) ^ C l b (R d ) for I < k - d/2. 

In this example L : H k (R d ) — s- H- k (R d ) is the operator 
L = (Id —A) k and the kernel K is the Green's function 
of L, 



K(x,y) = {2^)- d ' 2 J j—^\x- V \ k -2J k ^{\x-y\) 



Id 



J a (x) is the modified Bessel function of order a. Around 
x = the Bessel function behaves like J a (|a;|) ~ \x\ a 
and so 

\x~y\ k -i J d(\x-y\) ~ \x-y\ 2k ~ d around x-y = 0. 

2 

Thus K e C 2fe - !i - 1 ( K<ixKd ) K<ixd ); this wil1 be relevant 
in the case of landmarks. 

In the above example L was a scalar differential op- 
erator; it acted on each component of the vector field 
equally and was a multiple of the identity matrix. This 
is not always the case. For example the operator asso- 
ciated to the family of a-6-c-metrics is in general not 
scalar and the corresponding kernel is a dense (not 
sparse) matrix. 

In Sect. [7] the metric on Diff (R d ) was introduced 
by choosing a differential operator. Given an operator 
with appropriate properties, it is possible to reconstruct 
the space H. The reason for emphasizing the space T-L 
and the reproducing kernel is twofold: first, the induced 
metrics on Emb(M, R d ) and the space of landmarks 
have a simpler representation in terms of the kernel. 



Second, in the literature on LDDMM, e.g. in 129 , the 



starting point is the space T~L of vector fields and by 
presenting both approaches we show their similarities. 



8.3 The metric on Emb(M,R d ) 

Let G u be a right- invariant metric on Diff(]R d ). The 
induced metric on Emb(M, R d ) is defined via 

G*(h,h)= inf (X,X) H . 

Ao q—h 

To compute a more explicit expression for G^, we de- 
compose % into 

H™ Tt = {X : X\ q = 0} , Uf r = (H™ rt )^ . 



Then the induced metric is 
G' q i (h,h) = {X hoi ,X hoi ) n , 

with X £ X(R d ) any vector field such that Xoq = h. 
The horizontal projection does not depend on the choice 
of the lift, i.e., if X, Y g X(M d ) coincide along q, then 
x hor = yhc* We identify U\ OT with the G w -completion 
of the tangent space T q Emb(M, R d ). There are maps 

T q Emb(M, R d ) -> U 1 ™ -> C\ (Af, R d ) 



h 



i-> X 



hor ) 



X i-» Xoq 



The composition of these two maps is the canonical 
embedding T q Emb(Af, R d ) C£(M,M. d ). The space 
%k° r is again a reproducing kernel Hilbert space with 
the kernel given by 



K q : M x M 



ndxd 



, K q (x,y) = K(q(x),q(y)) . 



Thus we have identified the induced Riemannian metric 
G H on Emb{M,R d ) as 

G?(h,h) = (h,h) n nor. 



In this formula we identified "H^ or with vector fields on 
M with values in R d . 



8.4 Geodesic distance 

If the metric on Diff (R d ) is strong enough, then the in- 
duced Riemannian metric on Emb(M, R d ) has a point- 
separating geodesic distance function and we conjecture 
that the same is true for B e (M, R d ). 

Theorem 8.2 If the norm on % is stronger than the 
C^-norm. i.e 
we have 

ho - 00 < distE mb (go,gi) 



, H ^ C$, then for q Q ,q x € Emb(M,M d ) 



Proof For h G Emb(M,]R d ) and x G M let X € X(R d ) 
be any vector field with X o q = h. Then from 

\h(x)\ = \x( q (x))\<\\x\\ 00 <\\x\\ n , 

we see that ||/i||oo < Halloo and by taking the infimum 
over all X we obtain 

\\h\\l,<G%(h,h). 

Now for any path q(t) between go arid q\ we have 



qi(x)-q (x)= / d t q(t,x)dt, 
h 
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and thus 

\ qi (x) - q (x)\ < [ \d t q(t,x)dt\ < 
Jo 

Jo 



^G q(t) H(d t q(t),d t q{t j) dt = Len^ mb (g) . 



By taking the supremum over x G M and the infimum 
over all paths we obtain 



for the derivative of the kernel with respect to the first 
variable. The geodesic equation is given by 

d t q(t,x) =X(t,q(t,x)) 
8 t (p(t) ® vol 9 ) (t, x) = 

p(t, x) T K' q(t) (x, y)p(t, y) voF(y) ) ® vol 9 (a;) 



M 



Mo - gilloo 

as required 



< dist 



Emb 



(go, si) 



X(t, u) = / K(u, q(t, y))p(t, y) vol 9 (y) . 

JM 



See 86 for a derivation of the geodesic equation for 



For the geodesic distance on shape space we have a 
positive result for the space B e (S 1 ,M. 2 ) of plane curves 
and the family H — H k (R d ) of Sobolev spaces. A lower 
bound on dist^ e is given by the Frechet distance Q. 

Theorem 8.3 The geodesic distance on B e (S 1 ,M 2 ) of 
the outer metric induced by H = H k (R d ) with the op- 
erator L = (1 — AA) k for A > and k > 1 is bounded 
from below by the Frechet distance, i.e., for Qo,Q± G 
B e (S\]R 2 ) we have 

dist£7(<9o,<9i) < distf*(Q ,Qi)- 

Proof Take Q , Q x G B e (M, R d ) and let Q(t) be a path 
between them. Then by [86j Prop. 5.7] we can lift this 
path to a horizontal path q(t) on Emb(S' 1 , M 2 ). Then 



plane curves and 84 for the related geodesic equation 
on B e (M,R d ). 

The vector field X is not smooth but only X G H. 
Therefore it is not possible to horizontally lift geode- 
sies from Emb(M, R d ) to Diff(IR d ). One can however 
work in a suitable Sobolev completion of Diff (R d ). Then 
right-invariance of the Riemannian metric on Diff (R d ) 
implies the conservation of the momentum: 

ip(t)*p(t, .) is independent of t . 

From here we obtain via p(t, .) = ip(t)*p(0, .) that ip(t) 
satisfies the Lagrangian form of EPDiff, 



9t<p{t,x)= K(<p(t,x),<p(t,.)<p(t) m p(p,.). (18) 

JM 



See 56 for details on singular solutions of the EPDiff 



distl. (Q ,Qi)< lk(0) 



g(i)IU = 

= Len Emb (g) 



equations. Theorem |7.5| can be applied to show long- 



Lens (Q) time existence of solutions of ( 18 1 



and by taking the infimum over all paths we obtain the 
result. □ 

In order to generalize this result to B e (M, M. d ) one 
would need to be able to lift horizontal paths from 
B e (M,R d ) to Emb(M,M d ). A careful analysis of the 
induced metric G n in the spirit of [86] should provide 
such a result for a fairly general Sobolev-type metric. 



8.5 Geodesic equation 

The geodesic equation on Emb(Af, R d ) is most conve- 
niently written in Hamiltonian form in terms of the 
position q(t) and the momentum pit) = p(t) £§> vol 9 , 
where vol 9 = vol(g) = vol(g*<?). The momentum de- 
fines a vector-valued distribution with support in the 
image of q(t). The momentum p acts on X G X(E d ) by 

(Xoq(t),p(t))= f (X(q(t,x)),p(t,x))voV(x). 
M JM 

Let us introduce the notation 
K' q(je) (x,y)=D 1 K(q(t,x),q(t,y)) 



8.6 Curvature 

The representation of B e (M, R d ) as the quotient 
B e (M, R d ) = Diff (R d )/ Diff (K% 



was used in 84 together with an infinite dimensional 



version of O 'Neil's formula to compute an expression 
for the sectional curvature on B e (M, R d ). For details 
consult (84l Sect. 5]. 



9 The space of landmarks 

By choosing M to be the finite set M = {1, . . . , n} we 
obtain as Emb(M, R d ) the set of landmarks, i.e., the set 
of n distinct, labeled points in R d . Let us denote this 
space by 



C n (R d ) := {(q\ 



.,q n )\q k &R d ,q k ^q\k^j} 



Note that C n {R d ) is an open subset of R nd and thus it is 
the first example of a finite dimensional shape space in 
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this paper. As a consequence some of the questions dis- 
cussed for other shape spaces have a simple answer for 
the space of landmarks. The geodesic distance is guar- 
anteed to be point-separating, the geodesic equation is 
an ODE and therefore locally well-posed and due to 
Hopf-Rinow geodesic completeness implies metric com- 
pleteness. 

Remark 9. 1 We regard landmark space as the set of all 
labeled collections of n points in R d , i.e., the landmarks 
q = (q 1 , q 2 , . . . q n ), q = (q 2 , q 1 , . . . q n ) are regarded as 
different elements of C n (M. d ). One could also consider 
the space of unlabeled landmarks which would 

correspond to B e (M,Mr). It is sometimes called also 
configuration space. Since Diff (M) = S n is the symmet- 
ric group of n elements, we have = C n (R d )/ S n . 
The group S n is a finite group, therefore the projec- 
tion C n (M. d ) —> is a covering map and so for 
local properties of Riemannian geometry it is enough 
to study the space C n (R d ). 

Before we proceed we need to fix an ordering for 
the coordinates on M. nd . There are two canonical choices 
and we will follow the convention of [60] , A landmark 
q is a vector q — (q 1 , . . . , q n ) T £ C n (R d ) and each q l 
has d components q % — (q zl , . . . , q ld ) T . We concatenate 
these vectors as follows 



li 



,q ld ,a 2 \ 



2d 



nd\T 



(19) 



Riemannian metrics on C n (M. d ), that are induced 
by the action of the diffeomorphism group, have been 



studied in 60 78 83 and on the landmark space on 



the sphere in 151] . Other metrics on landmark space 
include Bookstein's thin-plate spline distance 



and Kendall's similitude invariant distance 64 . See 



18 



for an overview comparing the different approaches. 



9.1 The metric on C n {l 



As is Sect. 8.2 let the metric G H on Diff (R d ) be defined 



via a Hilbert space H of vector fields and let K be the 
reproducing kernel of H. As before Diff (R d ) can be any 
of the groups Diff c (R d ), Diff 5 (R d ) or Diff#= (R d ). The 
metric G H induces a Riemannian metric g-u on C n (R d ) 
and we can calculate it explicitly; see Thm |9.3| 

For the convenience of the reader we will repeat the 
definition of the distance function on £ n (M. d ) induced 
by the metric G H ; see Sect. [5] for the the more general 
situation of embeddings of an arbitrary manifold M in 
K d . Let E be the energy functional of the metric G H 
on the diffeomorphism group, i.e., 



E{v)= ThV)^ 
Jo 



dt. 



(20) 



The induced distance function of the action of the dif- 
feomorphism group on the landmark space is given by 



dist w (g,g)=inf{v/£H : ^( 9 l ) = , 



(21) 



where ip v is the flow of the vector field v and where 
the infimum is taken over all sufficiently smooth vector 
fields v : [0, 1] — > X(M d ). Given a solution v of the above 
minimization problem, the landmark trajectories q l (t) 
are then given as the solutions of the ODE 



q i (t)=v(t,q i (t)), 



1. 



We will now define a Riemannian metric on the fi- 
nite dimensional space £ n (M. d ) directly and we will see 
that it is in fact induced by the metric G H on the dif- 
feomorphism group. For a landmark q we define the 
matrix 



9h(i) = 



where K 



(K^.q 1 ) ■■■ K{q\q n y 



XK^.q 1 )-- - K(q n ,q n ) / 



vndxnd 



(22) 



xl^ R dxd is the kernel of H. That g H 
defines a Riemannian metric on C n (W l ) can be easily 
shown using the properties of the kernel K. 

The metric gn defines, in the usual way, an energy 
functional directly on the space of landmark trajecto- 
ries, 

E(q(t))= f q(t) T g n (q(t))q(t)dt, 
Jo 



(23) 



and one can also define the induced distance function 
of g as 



?(*) 



dist H (g, q) = inf V E{q{t)) : g(0) = q, q(l) = q 



(24) 



where the infimum is taken over all sufficiently smooth 
paths in landmark space q : [0, 1] — > C n (M. d ). 



It is shown in 83 Prop. 2] that the minimization 



problems ( 20 ) and ( 23 ) are equivalent and that the in- 



duced distance functions are equal: 

Theorem 9.2 (Prop. 2, |83|) Let v be a minimizer 
of the energy functional ( |20[ ). Then the trajectory q{t), 
which is obtained as the solution of the system of ODE's 

q'i^^v^q^t)), i=l...,n, 



minimizes the energy functional (23) and E(v) — E(q). 



On the other hand, if q{t) is a minimizer of the energy 



functional ( 23 I define the vector field 

n 

v{t,x)=Y J Pi(t)K{x,q i {t)) 

i=l 



(25) 
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with the momenta pi : [0, 1] 
p i (t)=K(q(t))- 1 q i (t). 



given by 



(26) 



Then the vector field v is a minimizer of the energy (|20j) 
and we have E(q) = E(v). 

Thus we have: 

Theorem 9.3 IfH then the induced metric gu 

on C n (M. d ) is given by 

(K(q\q x ) ■■■ K(q\q n Y 



9h{q) = 



r>ndxnd 



I*,? 1 )-' K(q n ,q n ), 



(27) 



where K G C k {R d x M d ,M dxd ) is the kernel ofU. We 
have g H € C k (W ld ,R ndxnd ). 

We will discuss the solutions of the minimization prob- 



lem (231 in Sect. 9.3 



Remark 9.4 Note that in the articles 82 83 the co- 
ordinates were ordered in a different way. Given q = 
(q 1 , . . . , q n ) they flatten it as 



(I 1 



,q na ) 



d\T 



If the kernel K(x,y) of H is a multiple of the iden- 
tity matrix, i.e., K(x,y) = k(x, y) Id^xd for a scalar 
function k, then the matrix g<n(q) is sparse and these 
coordinates allow us see the sparsity in an elegant way, 



9n(l) = 



(k{q) 
k(q) 



° \ 





V ••• k(q)J 
Here k(q) denotes the n x n-matrix (fc(<A <7"'))i<ij<n- 

9.2 The geodesic equation 

The geodesic equation can be deduced from the equa- 



tion in the general case Emb(M, M. d ); see Sect. 8.5 

Theorem 9.5 IfH <-> C\, then the Hamiltonian form 
of the geodesic equation of the metric gu on C n (M. d ) is 
given by 



N 
3=1 

N 

Pi = - J £pJ(diK)(q i ,qj)p j 



(28) 



with Pi(t) = K(q(t)) q l (t) the vector valued momen- 
tum. 



For scalar kernels this system has been studied in the 



articles [78j|S3j ; see also the PhD-thesis of Micheli [82 
Two examples of a two-particle interaction can be seen 
in Fig(8] 

Remark 9.6 A different possibility to derive the above 
geodesic equation is to consider directly the Hamil- 
tonian function of the finite dimensional Riemannian 



manifold (£ n (R d ) , g n ) . Following [61] Eqn. 1.6.6] it is 
given by 



Ham(p, q) = -p T g{q) l p 



N 



Then the geodesic equations (28) are just Hamilton's 
equation for Ham: 



d Ham 

dpi 



Pi 



d Ham 







m 


Willi 'JU 







Fig. 8 Two geodesies in C 2 (R d ). The grid represents the cor- 
responding diffeomorphisms. On the left-hand side both land- 
marks travel in the same direction, and the two paths tend 
to attract each other. On the right hand side they travel in 
opposite directions and thus try to avoid each other. Original 
image published in [83] . 



Remark 9. 7 We can regard a geodesic curve of land- 
marks as a soliton-like solution of the geodesic equa- 
tion on Diff(IR <i ) where the corresponding momentum 
is a linaer combination of vector valued delta distribu- 
tions and travels as such. 



9.3 Completeness 

As a consequence of the global well-posedness theorem 

we can 



on the full diffeomorphism group - Thm. 7.5 
deduce the long-time existence of geodesies on Land- 
mark space. To do so we solve the geodesic equation 



( 18 ) on the diffeomorphism group for a singular initial 
momentum p(0, x) = Pj°~( x — Q J )- Then the land- 

mark trajectories are given by q l (t) — ip{t, g I (0)), where 



Lp e Diff(Ifr) is the solution of (18). 
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Theorem 9.8 IfH C\, then the Riemannian man- 
ifold LC" (R d ), gu) is geodesically complete. 

A consequence of this theorem is that two landmarks 
will never collide along a geodesic path. For finite di- 
mensional Riemannian manifolds with a metric that is 
at least C 2 the theorem of Hopf-Rinow asserts that the 
notions of geodesic completeness and metric complete- 
ness are equivalent. 

Corollary 9.9 If U ^ C 2 , then (£ n (R d ), dist w ) is a 
complete metric space. 

For a C 2 -metric g-u one can use once more the theo- 
rem of Hopf-Rinow to show the well-posedness of the 
geodesic boundary value problem. 

Corollary 9.10 (Prop. 1, (83)) IfH <-> C h 2 then for 
each pair of landmarks q,q £ C n (R d ) there exists a min- 
imizer q(t) € C 1 ([0, 1], of the energy functional 

(23}. 




Fig. 9 Sectional curvature of £ 2 (R),gx as a function of the 
distance between the landmarks qo,qi- Here K was the Gaus- 



sian kernel K(x) = exp(— |). Original image published in 83 



In fact the existence of minimizers to the boundary 
value problem on landmark space can be proven un- 
der even weaker smoothness conditions on the metric 



g-u, see 60 Sect. C] 



9.4 Curvature 



We see from ( 22 ) that the expression for the co-metric 
g^, 1 is much simpler than that for g-u- In the article 83 
the authors took this observation as a motivation to 
derive a formula for the sectional curvature in terms 



of the co-metric, now called Marios formula; see 83 
Thm. 3.2]. Using this formula they were able to calcu- 
late the the sectional curvature of the landmark space 
(£ N (R d ),g n ); see [831 Thm. 9]. We will not present 



these formulas in the general case but only for the spe- 
cial case of two landmarks in M: 

Theorem 9.11 (Prop. 23, |83|) The sectional cur- 
vature onC 2 (M) depends only on the distance p = \q—q\ 
between the two landmarks q,q. For a metric g-u, with 
reproducing kernel K £ C 2 (M. d , M. d ), it is given by 



K(0) - K(p) 
HP) = K(0) + K( P ) K {P) 



2K(0)-K(p) a 
(^(0)+^( P )) 2 (> 



10 Universal Teichmuller space as shape space 

Here we sketch how Diff(S' 1 )/PS , L(2, R) parametrizes 
the shape space of simple closed smooth plane curves 
modulo translations and scalings and discuss the associ- 
ated Riemannian metric, called the Weil-Peterson met- 
ric. This metric has nonpositive curvature, is geodesi- 
cally complete, and any two shapes can be connected 
by a unique minimal geodesic. There exist soliton-likc 
solutions which are called teichons which are given by 
finite dimensional Hamiltonian system. They relate to 
geodesies of shapes like landmarks do to geodesies in 
Diff (R d ); see Sect. 9.7 This theory and the correspond- 



ing numerical analysis has been developed in 111 [| 11 2 
The use of teichons has been developed in |74| . 

Given a 1-dimensional compact smooth and con- 
nected submanifold T in M 2 = C inside the Riemann 
sphere C = C U {oo}, we consider its interior T int and 
its exterior T cxt which contains oo; these are smooth 
2-manifolds with boundary. Let Di nt and D ext denote 
the unit disk and the exterior of the unit disk respec- 
tively. By the Riemann mapping theorem there exists 
a smooth conformal map ^j nt : D; nt — ¥ T; n t, unique up 
to replacing it by ^j„t o A for a Mobius transformation 



A{z) 



bz 



with 



£ PSU(1,1) PSL{2, 



For a Gaussian kernel K a plot of the curvature de- 
pending on the distance between the landmarks can be 
seen in Fig|9j 



Likewise we have a conformal map between the exte- 
riors #ext : D ext — > Text which is unique by the re- 
quirement that <?ext(oo) = oo and ^ xt (oo) > 0. The 
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resulting diffeomorphism 



: S 1 -> S 1 



projects to a unique element of DiS + (S 1 )/PSL(2,R) 
(here S 1 is viewed as P 1 (M)). It is called the fingerprint 
of r. Any coset \P.PSL(2,M.) comes from a shape r, 
and two shapes give the same coset if they differ by a 
Mobius transformation in Aut(C) which fixes oo and 
has positive derivative at oo; i.e., by translations and 
scalings. 

One can reconstruct the shape r from the finger- 
print W.PSL(2, R) by welding: Construct a Riemann 
surface by welding the boundaries of D; nt and D cxt via 
the mapping W. The result is confomally equivalent to 
the Riemann sphere and we use a conformal mapping 
<P from the welded surface to the sphere which takes oo 
to oo and has positive derivative at oo. Then _T equals 
<? _1 (5 1 ) up to a translation and a scaling of C. An 
efficient numerical procedure for welding is described 
in (mJ[lT2 . 

The quotient T := Diff {S 1 )/PSL{2, R), also known 
as universal Teichmiiller space, is naturally a coadjoint 



orbit of the Virasoro group (see Sect. 7.2) and as such 
it carries a natural invariant Kahler structure; see 69 



The corresponding Riemann metric can be described 
as follows. For u e X(S 1 ) = C°°(S X ) we consider the 
Fourier series u{9) = X)nez a " e with a n = a_„ and 
the seminorm 



2 

WP 



\n 3 - n\\a n Y 



The kernel of this seminorm consists of vector fields of 
the form aie-' ie + a + aie ie ; i.e., ker(||-|| WP ) = $1(2, R). 
So this gives an inner product on the tangent space at 
the base point of T. This norm can also be defined by 
the elliptic pseudifferential operator L = H(dg+dg) via 
IM| = J s i L(u).ud9, where the periodic Hilbert trans- 
form T-L is given by convolution with cotan(|). The 
inverse of L is convolution with the Green's function 



K{9) = £ 



Arid 



|n|>l 



(I - cos 9) log (2(1 - cos 0)) + - cos - 1. 



According to Sect. 7.2 <p(t) G Diff(5' 1 ) projects to a 
geodesic in T if and only if the right logarithmic deriva- 
tive u{t) = d t <p(t) o (pit) -1 satisfies 

L{u t ) = — ad* (Lit) or 

{Lu) t + u.{Lu)g + 2ug.(Lu) = 



and u(0) has vanishing Fourier coefficients of order — 1, 
0, 1. We call m — Lu £ (X(S 1 )/sl(2, R))' the momen- 
tum, with u = G * m. The Weil-Petersen metric de- 
scribed by L is a Sobolev metric of order 3/2. The ex- 
tension to the corresponding Sobolev completions has 
been worked out by [49) . 

If we look for the geodesic evolution of a momentum 
of the form 



N 



^2pjS{9-qj), so that v = ^2,pjG(e -qj), 



3=1 



is a finite combination of delta distributions, which lies 
outside of the image of L : X(S 1 )/sl — > (X(S 1 )/sl)' , we 
see that the evolution of the parameters qj,Pj is given 
by the Hamiltonian system 



Pk 



PkTrLiPjG'iQk - Qj) 



These solutions are called Teichons, and they can be 
used to approximate smooth geodesies of shapes in a 
very efficient way which mimics the evolution of land- 
marks. The disadvantage is, that near concave parts of 
a shape the teichons crowd up exponentially. An exam- 
ple of such a geodesic can be seen in Fig 



10 



sec 



74 



and 75 for more details. 




Fig. 10 Evolution of an 8-Teichon from the circle to a 
Donald-Duck-like shape. Positions of individual 1-Teichons 
are marked by asterisks. Original image published in [74] . 



11 The space of Riemannian metrics 

Let M be a compact manifold without boundary and 
dim(M) = rn. In this part we describe the Riemannian 
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geometry on Met(M), the manifold of all Riemannian 
metrics on M. The L 2 -metric on Met (AT) is given by 



G E (h,k) = [ Tr (g-'hg-'k) vol(g) , 

JM 



with g e Met (A/) and h,keT g Met(AT). Each tangent 
vector h is a bilinear form h : TM x M TM — > E, that 
is interpreted as a map TM — > T*M. This metric has 



been introduced in 41 and is also known as the Ebin- 



metric. Its geodesic equation and curvature have been 
calculated in 47,50 , and the induced distance function 



and metric completion have been studied by Clarke 27 
30l 



Similar to Riemannian metrics on immersions, So- 
bolev metrics of higher order and almost local metrics 
can be defined using a (pseudo differential) operator L 
acting on the tangent space of Met(M). To be more 
precise, let 

L : TMet(Af) -» TMet(M) 

be a smooth base-point preserving bundle isomorphism, 
such that for every g £ Met(M) the map 

L g : T g Met (A/) ->■ T g Met(M) 

is a pseudo differential operator, that is symmetric and 
positive with respect to the metric G E . Then we can 
define the metric G L by 

G L (h, k)= [ Tr (g^Lgih^k) vo\(g) . 

JM 

Let us also assume, that the operator L is invariant 
under the action of Diff(M), i.e., 

ip*{L g h) = L v * g (ip*h) . 

Then the metric G L induces a Riemannian metric on 
Met(M)/Diff (M) where Diff (M) denotes the group 
of all diffcomorphisms that are homotopic to the iden- 
tity. In relativity theory the semi Riemannian analogon 
of Met(M)/ Diffo(Af) is called super space, since it is 
the true phase space of Einstein's equation. 
An example for an operator L is 



L g h = h+(A 9 yh, 



l>0. 



The resulting metric G L , which is a the Sobolev metric 
of order I, has been introduced in [12]. Other metrics, 
that have been studied include conformal transforma- 



tions of the L 2 -metric 12 31 



G*(h, k) = #(Vol 9 ) f Tr (g-'hg^k) vol(g) 

JM 



with <P € C°°(IR>o,R>o) and scalar curvature weighted 
metrics [12] , 

G*(h, k) = f <T>(Scal 9 ) Tr (g^hg^k) vol(g) , 

JM 

with <2> G C°°(]R,]R > o). 

The main focus of the section will be on the L 2 - 
metric. 



11.1 Connections to Teichmuller theory and 
information geometry 

Our main motivation to consider the space of all Rie- 
mannian metrics in this article lies in it's possible ap- 
plication to shape analysis of surfaces as explained in 
Sect. [TTJ see also [59 



Another motivation for the study of the L 2 -metric 
on the manifold of metrics can be found in its con- 
nections to Teichmuller theory. Let M be a Ricmann 
surface of genus greater than one. Then the L 2 -metric, 
restricted to the space Meti of hyperbolic metrics, in- 
duces the Weil-Peterson metric on Teichmuller space 
Meti(M)/Diff (M). This is described for example by 



Fischer and Tromba 1461 or Yamada 124 125 



A surprising connection can be also found with the 
field of information geometry, since the L 2 -metric de- 
scends to the Fisher-Rao metric on the space of volume 
densities. To understand this connection we will con- 
sider the Riemannian metric on Diff(A/) induced by 
G E . For a fixed metric g G Met(M) we introduce the 
map: 

Pull : Diff(M) Met(M), tp ^ ip*g . 

Now we can define a metric G Pu11 on Diff(AT) as the 
pullback of the L 2 -metric under the map Pull, i.e., 

G? ul \h, k) = G E {T V Pull h, T v Pull fc) . 

This mapping and the induced metric on Diff(M) for 
a variety of metrics on Met (AT) is studied in [jo] . The 
metric G Pu11 is invariant under the left action by the 
group Diff A1 (M) of volume-preserving diffcomorphisms 
and the metric induced on the quotient space 

Dens(A/) = Diff M (Af)\ Diff (M) 



of densities is the Fisher-Rao metric; see 97 Thm. 4.9]. 

Another possibility, to see the connection to infor- 
mation geometry was implicitly presented in [32] . There 
the authors consider the subspace of Kahler metrics in 
a fixed Kahler class - assuming that M admits a Kahler 
structure. Then the Ebin metric induces the so-called 
Calabi geometry on the space of Kahler metrics. It was 



then observed in 67 that this space is, via the Calabi- 



Yau map, isometric to the space of volume densities 
with the Fisher-Rao metric. 
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11.2 Geodesic distance 

In contrast to to the spaces of immersions, submani- 
folds and the diffeomorphism group, the L 2 -metric on 
Met(M) induces a point-separating geodesic distance 
function. 

Theorem 11.1 (Thm. 18, (28)) The L 2 -metric in- 
duces a point- separating geodesic distance function on 
Met(M). 

Remark 11.2 Note that this result also holds for all 
metrics, that are stronger than the L 2 -metric, i.e., 

Gf(h,h)<CG g (h,h), 

with a constant C > 0, independent of g. This applies 
in particular to almost local metrics, if the function 
is bounded from below by <P > C > 0, as well as to 
most Sobolev-type metrics. 

Fix a Riemannian metric g on M. For each x G M 
denote by Met(M) a = S\T*M the space of symmet- 
ric, positive definite (2) -tensors at x. Then for b, c € 
T a Met(M) x the expression 

lx,a{b,c) = Tr(a~ 1 ba~ 1 c) ^dctg(x)~ 1 a 

defines a Riemannian metric on the finite-dimensional 
manifold Met(M) x . Denote by d x the induced geodesic 
distance function and define the following distance on 
Met(M), 



^2(90,91) 



d x {go(x),gi(x)) 2 vol(g) . 



The following theorem states that computing the geode- 
sic distance on Met(M) with respect to the L 2 -distance, 
is equivalent to summing point-wise geodesic distances 
on Met(M) x . 

Theorem 11.3 (Thm. 3.8, [3l]) The geodesic dis- 
tance induced by the L 2 metric and the distance Q2 
coincide, i.e., for all go,gi € Met(M), 

dist B (.g ,5i) = tt 2 (g ,gi) ■ 

Similar as in the case of the G^-metric and the So- 
bolev metrics on the space of immersions the square 
root of the volume is again a Lipschitz continuous func- 
tion. 

Theorem 11.4 (Lem. 12, |28|) The induced geodesic 
distance induced by the L 2 -metric satisfies the inequal- 
ity 



for any measurable set F C M. Here Yol(F,g) denotes 
the volume of F C M with respect to the metric g. 
This implies the Lipschitz continuity of the map 



Wol: (Met(M),disC ot(M) J -> R> . 

On the other hand we also have the following upper 
bound for the geodesic distance. 



Theorem 11.5 (Prop. 4.1, [30]) For the L 2 -metric 
the geodesic distance is bounded from above by 



dist £ ( 50 ,ffi) < C(m) (yVol(F, 3o ) + y/Vol(F, 9l )) , 
where E denotes the support of g\ — go 



F = {x e M\g {x) ^ 9l {x)} , 

and C(m) is a constant depending only on the dimen- 
sion of M . 

The above corollary implies that the set Met M (M) 
of all Riemannian metrics having a total volume less or 
equal than fi has a finite diameter with respect to the 
L 2 -metric. 



11.3 The geodesic equation 

The Christoffel symbols for the L 2 -metric were first cal- 
Sect. 4]. Subsequently Freed and Groisser 



culated in 



41 



47 and Michor and Gil-Medrano |50] computed the 
geodesic equation and found explicit solution formu- 
las. The geodesic equation for higher order Sobolev 
type metrics and Scalar curvature metrics can be found 



in [12] and for volume weighted metrics in 12 32 



The geodesic equation for the L 2 -metric decouples 
the time and spatial variables, i.e., instead of being a 
PDE in (t, x), it is only an ODE in t. 



Lemma 11.6 (Sect. 4, |41|) The geodesic equation 
for the L 2 -metric is given by the ordinary differential 
equation: 

9tt = \ Tr(g~ 1 gtg~ 1 gt)g + gtg~ x gt - ^ Tr( 5 _1 5t ) 3t . 

There exists an explicit solution formula for this ODE. 



Theorem 11.7 The geodesic starting at go € Met(M) 
in the direction ofh€ T go Met(M) is given by the curve 



VVol(F,<7o)-VVol(F, ffl ) 



< 



dist 



Mct(M)(So,5i) g(t) = ffoe «(Oid+6(t)Ho , 
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where Hq is the traceless part of H := g h, i.e., Hq = 
H - '^p-Id, and where a(t) and b(t) G C°°(M) are 
defined by 

a{t) = I log ((1 + f Tr(ff)) 2 + § Tr(ff 2 )t 2 ) 

4 arctanf ^S^ ),IV(j?g)^0 

Tr(i/„ 2 ) = 0. 



1+f Tr(ff)' 



Here arctan is taken to have values in (— -j, 5) for the 
points of the manifold where Tr(H) > 0, and on a point 
where Tt(H) < we define 



arctan 



^mTr(iJ 2 )A 
4 + tTr(F) J 

arctan in [0, f ) for t G [0, — Bpn ) 

fort = - xFpry 
arctan in (f,7r) /or t 6 (- Tr(g) ' °°)- 



Let TV' 1 := {a; G M : H (x) = 0}, and if N h £ Zei 
t h := \ni{Tv(H){x) : x G iV 1 }. TTien t/ie geodesic g(t) 
is defined for t G [0, oo) i/ A' 1 = or i/ i ft > 0, and it 
is only defined for t £ {0, — jk) if t h < . 

These formulas have been independently derived by 



Freed and Groisser 47 and Michor and Gil-Medrano 
50 . A similar result is also available for the metric 
with #(Vol) 



l . 

Vol' 
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Remark 11.8 The geodesic equation for higher order 
metrics will generally not be an ODE anymore and ex- 
plicit solution formulas do not exist. Nevertheless, it 
has been shown that the geodesic equations are (lo- 
cally) well-posed, assuming certain conditions on the 



operator L defining the metric; see 12 . These condi 



tions are satisfied by the class of Sobolev type metrics 
and conformal metrics but not by the scalar curvature 
weighted metrics. 



11.4 Conserved quantities 

Noether's theorem associates to any metric on Met(M), 
that is invariant under pull-backs by the diffcomor- 
phism group Diff(M), for each X G X(M) the quantity 

G g (g t ,Cx(g)) = const, 

which is conserved along each geodesic g(t). Here (x 
is the fundamental vector field of the right action of 
Diff(M), 

Cx(g) = C x g = 2SymV g (g(X)), 



and Sym \7 9 (g(X)) is the symmetrization of the bilinear 
form (Y, Z) h> \7 9 Y g(X, Z), i.e., 

SymV° (g(X))(Y,Z) = \ (V^g(X, Z) + V 9 z g(X, Y)) . 

If G g (gt,Cx(g)) vanishes for all vector fields X G X(M) 
along a geodesic g(t), then g(t) intersects each Diff (M)- 
orbit orthogonally. 



11.5 Completeness 

The L 2 -mctric on Met(M) is incomplete, both metri- 
cally and geodesically. The metric completion of it has 
been studied by Clarke in 27 30 31 . To describe the 



completion let Met/ denote the set of measurable sec- 
tions of the bundle S> T*M of symmetric, positive 
semi-definite ( 2 )-tensors, which have finite total vol- 
ume. Define an equivalence relation on Met/ by iden- 
tifying go ~ <?i, if the following statement holds almost 
surely: 

9o(x) 7^ 9i{x) =^ both <7i(x) are not positive definite. 

In other words, let D = {x : go(x) ^ gi(x)} and Ai — 
{x : gi(x) not pos. dcf.}. Then 

go ~ gi & D \ (Ao n A\) has measure . 

Note that the map Met(M) ^ Met/ / ~ is injective. 

Theorem 11.9 (Thm. 5.17, [27]) The metric com- 
pletion of the space (Met(M),dist fi ) can be naturally 
identified with Met/ /~. 

In the subsequent article |31| it is shown that the 
metric completion is a non-positively curved space in 
the sense of Alexandrov. 



Theorem 11.10 (Thm 5.6, [31]) The metric com- 
pletion Met(Af) o/Met(M) with respect to the dist B - 
metric is a CAT(0) space, i.e., 

1. there exists a length-minimizing path (geodesic) be- 
tween any two points in Met(Af) and 

2. (Met(A-f),dist B ) is a non-positively curved space in 
the sense of Alexandrov. 



11.6 Curvature 

For the L 2 -metric, there exists a comparably simple 
expression for the curvature tensor. 
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Theorem 11.11 (Prop. 2.6, 1 50] ) The Riemannian 
curvature for the L 2 -metric on the manifold Met(M) 
of all Riemannian metrics is given by 

g- x R g (h,k)l=\[[H,K],L] 



16 
1 

16 
1 



(Tr(KL)H - Tr(HL)K) 

(Tr(£f) Tr(L)K - Tr(K) Tr(L)H) 



+ — (Tr(K) Tr(HL) - Tr(H) Tt(KL)) Id . 
16 

where H = g~ 1 h, K = g k and L = g~ x l. 

In the article [47] the authors have determined the 
sign of the sectional curvature: 

Theorem 11.12 (Cor. 1.17, (47)) The sectional cur- 
vature for the L 2 -metric on the manifold Met(Af) of 
all Riemannian metrics is non-positive. For the plane 
P{h, k) spanned by orthonormal h, k it is 



.Met 



(P(h, k)) = J ^ ( Tr(HK) 2 - Tr(H) 2 Tr(K) 2 
+ 1 -Tr((g- 1 [h,k}) 2 )vol(g), 



where H — g h and K = g 1 k. 



In 32 it is proven that this negative curvature carries 
over to the metric-completion of (Met(M), G E ), as it is 



a CAT(O) space; see Lem. 11.10 
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